1176 00189 7074 ^ 




TECHNICAL MEMORANDUMS 
NATIONAL ADVISORY COMMITTEE EOR AERONAUTICS 




No. 1026 



THERMAL THEORY. OE COMBUSTION AND EXPLOSION 
III. ..THEORY OP NORMAL. ELAME PROPAGATION 
By N. N. Semenov 

Progress of Physical Sc ience: (U. 5 . S .R . ) 
Vol. 24, No. 4, 1940 



|f A € A UHRAKT 

M$$$$f MWffiOWM' AERONAUTISAk- 
LABQRATOR¥ 
feaagtey Field, Vsu .^mS 

if* 



Washingt on . 
September 1942 



HAT IOH A L ADVISORY COMMITTEE EOR AERONAUTICS 



TECHNICAL MEMORANDUM NO. 1026 



THERMAL THEORY OE COMBUSTION AND EXPLOSION 
III. .THEORY OE NORMAL ELAME PROPAGATION* 
By N. N. Semenov 

1. EXPERIMENTAL DATA ON ELAME PROPAGATION 



On igniting a combust ible gas mixture at any point by 
means of a spark or hot body the process of combustion will- 
be propagated throughout the mass of the mixture. This is 
what occurs in the explosion of methane in coal mines and 
of the fuel mixture in the internal combustion engine. If 
this phenomenon is observed through a glass tube filled 
with the combustible mixture and ignited by a spark at one 
end, it will readily be found that along the tube a narrow 
flame front is propagated, which separates the already 
burning gas (behind the front) from the fresh as yet un- 
ignited gas (ahead of the front). 

In the above cases the flame moves through an initially- 
stationary gas. In many cases, however, the process of 
combustion takes place in a combustible gas which is moving 
toward the ignition source. A stationary flame will then 
be obtained such as that in the usual Bunsen burner. In 
each case we are dealing with the same phenomenon of the 
motion of the flame front with respect- to the unburned 
gas. In the combustion of a • homogeneous mixture the com- 
bustion proceeds with groat intensity because the velocity 
of the reaction at the flame temperature is very high. The 
process occurs with much less intensity when the combusti- 
ble gas is not completely mixed with air, as is the case, 
for example, in the burning of a naphtha torch flame issuing 
from a nozzle or in the combustion of a gas burner without 
air opening. In these cases the rate of combustion Is 
limited not by the speed of reaction but by the rate 'of 
the mixing through diffusion of the molecules of the com- 
bustible and the oxygen. In Diesel engines and forced 
furnaces it is attempted to accelerate this process of 
mixing of the reacting gases by constructional means that 
will insure the energetic turbulence of the gas. Not- 
withstanding various mathematical difficulties the theory 



Progress of Physical Science." (u.S.S.R.), vol. 24, no. 4, 
1940, pp. 433-486, 



2 _> NACA Technical Memorandum No. 1026- 



of purely diffusive flames is fundamentally clear since 
the rate of combustion is determined by simple diffusive 
processes. Much less understood is the fundamental proc- 
ess of combustion of the combustible gas mixtures which 
was mentioned earlier. In the present paper \*e shall con- 
cern ourselves precisely with this problem to which our 
discussion will be limited. We may note, incidentally, 
in a number of cases the combustible fuel gas is always 
"mixed". This is the case when the flame is associated 
with the decomposition of any single given endothermic 
substance and proceeds without oxygen or air (decomposi- 
tion of C1 3 0, ozone, vapors of explosive gases, etc.). 
By the term combustion we shall denote the propagation 
of the flame in the process of combining with oxygen or 
during any other chemical processes where a certain 
amount of energy is liberated and a flame may exist. 

The temperature of the gases immediately behind the 
flame front can be readily computed from the thermodynamic 
data if the heat losses in the flame zone are not large. 
This temperature for various compositions of the mixture 
fluctuates generally between 1500° and 3000° K. Often (oh 
account of dissociation) complete combustion does not take 
place in the flame front and a part of the substance may 
complete its combustion behind the flame front. The Bun- 
sen burner generally operates with excess fuel mixtures. 
The most intense combustion takes place in the inner cone 
which forms a thin flame cap. In this zone the fuel burns 
to CO, C0 3 , H 3 , and H 2 0, the relative concentration of 
which is determined by the equilibrium of the hydrogen 
gas at the combustion temperature. The products CO 
and H 3 are burned by the diffusion of the oxygen of 
the air within the wide outer cone of the burner. In 
this zone the process is determined not by the speed of 
reaction but by the diffusion. 

The two processes in the Bunsen burner may be sepa- 
rated, for which purpose the burner is surrounded by a 
glass tube held in place by a plug, as shown in figure 1. 
By this means the cones can be separated. The inner cone 
(a) as before forms a kind of thin flame cap over the 
end of the burner. Since the glass tube becomes filled, 
with the products of combustion intermingled \irith the 
products of incomplete combustion and dissociation of the 
fuel, there is no access of the oxygen of the air to the in 
side cone and therefore no outer cone of the flame is here 
formed. At the exit to the air, however, these products 
can burn by mixing with the oxygen and thus the outer 
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cone (Id) is displaced at the mouth of the tube and 
the two cones separated. 

¥e shall thus concern ourselves with the process of 
combustion in homogeneous mixtures taking place in a very- 
thin zone of the flame front. (Later we shall see that the 
thickness of this zone is of the order of 10 ® era, ) This 
thin zone is sufficient, however, for the reaction proc- 
esses to take place, the heat given out being expended 
in heating the fresh cool gas ahead of the flame front 
to a high temperature. Due to the fact that the gas on 
one Bide of the zone consists of the products of reaction 
and on the other side of the initial substances energetic 
diffusion processes will take place in this zone. These 
in sum are the phenomena that occur in the flame front. 
Our problem is, on the ba.sis of a study of these phenomena, 
to obtain quantitative expressions for the speed of com- 
bustion of the substances in the flame front. 

Due to the fact that the zone within which the com*- 
bustion takes place is very thin, external, artificially, 
produced turbulence and even the motion of the gas arising, 
from the combustion process, will lead only to a curvature 
of the surface of the flame front and an increase in its . 
area but will not disturb the structure of the zone itself, 
(They will not, for example, complicate the processes of 
diffusion and heat transfer in the zone,) Thus any kind 
of motion or turbulence in the gas only increases the 
area of the flame front but does not change the' charac- 
teristics of any element of the front area. Particularly, 
on a unit area of the flame front of a given mixture the 
same quantity of fuel per second v m will always be burned. 
Therefore the total quantity of fuel burned per second over 
the entire flame front area S will be v m S. Any curvature 
of the flame front due to gas flows or turbulence will 
therefore increase the total quantity of burned fuel in 
proportion to the increase in the flame front area S. 
There thus enters in the theory of combustion for each 
given mixture a certain constant , namely, the mass rate of 
combustion v m equal to the number of grams of mixture 
burned per second per unit area of the flame front. 

In burning a definite quantity of fuel per second, - 
the flame front is displaced with respect to the initial 
gas as fresh quantities of the latter are consumed by the 
flame. It is not difficult to establish a relation between 
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the mass rate of combustion and the linear velocity of 
the displacement of the flame front. In order to secure 
the combustion of v m dS grams of the mixture it is 
necessary that this quantity of the mixture be brought 
up to -the - element dS of the flame front either by 
displacement of the flame front or by having fresh gas 
supplied to the stationary flame front. For this pur- 
pose the linear velocity of displacement of a given 
element of the flame front with respect to the unburned 
gas in the direction normal to the flame front surface 
at a given position should be equal to v 0 = v m /p 0 where 
p is the density of the initial cold gas. This linear 
velocity v Q is denoted as the "normal velocity of flame 
propagation. "* It is a characteristic constant for the 
combustion process of every mixture and does not depend 
on the hydrodynamic conditions under which the combustion 
occurs . 

If the flame front in a tube were plane (or a sphere 
normal to the axis of the tube) the displacement velocity 
of the flame front (with ignition at the open end of the 
tube) would exactly be equal to the normal speed of prop- 
agation. Actually the burned gases on expanding flow 
out through the open end of the tube into the atmosphere 
according to the law of laminar flow (that is, with veloc- 
ities varying with the distance from the wall). At the 
wall itself the flow velocity of the burned gases is equal 
to zero. It can readily be shown that the burning gases 
at the flame front edges should flow at an angle to the 
latter. A detailed though qualitative analysis of the 
hydrodynamics of the flow of the gases in this case (see 
Jost reference 1) leads to the conclusion that the plane 
front is unstable and that the front should assume a 
curvature. The stable form will be th?.t of a cap which 
bulges in the direction of motion of the flame. On figure 
2 is shown a number of instantaneous flame-front photo- 
graphs separated by equal time intervals, the front being 
propagated in a vertical tube (d = 5 cm, mixture of 58 
percent CO and 42 percent air, pictures by Barsky). 
For flame propagation in a horizontal tube (if the latter 
is not too narrow) there will also be a convexity due 
to gravity which makes the front unsymmetrical with 
respect to t he tube axis (fig* 3) . As long as there 
is no change in the hydrodynamic conditions the shape of the 
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flame front will remain 'constant, hence also the quantity 
v m S of '.gas burned per second and also the velocity of 
displacement of the flame front as a whole* This veloc- 
ity v 0 i) S will^evidently be equal to v m S/TTr s p 0 * where 
r is the radius' 1 of the tube and S as before the total 
'flame-front area, whence v 9tg = v Q S/tt r 3 , that is, the 
ratio of the observed velocity of propagation in a cer- 
tain tube to the normal velocity is as the area of the 
flame front to that of the tube. Tests show, in fact, that 
for a distance .0.5 - .1 meter (depending on the total length 
of the tube) the velocity of propagation of the flame in the 
tube 'remains constant. This is seen from figures 2, 3, and 
4, which shows a photograph of the flame propagation obtained 
on a rotating film. (If w is the linear speed of motion 
of the film tan" cp = v 0 -^ s /w.) - The constancy of the angle 
cp, that is, the* rectilinear trace of the flame indicates 
the constancy of the propagation velocity. By both of 
these photography methods the value of the observed veloc- 
ity of propagation can be obtained. 

TABLE I 



Velocity of flame propagation in 10- percent mixtures of 
methane with air in tubes of various diameters and com- 
putation of the normal velocity of propagation 



Diameter 


Direction 


Observed velocity 


Area of 


normal 


of tube 


of flame 


of propagation 


■ flame 


veloc ity 


(cm) 


propagation 




front 






( cm/ sec) 


( cm 3 ) 


(cm/ sec) ' 


10 


Horizontal 


111 


300 


29 


10 


Do 


71 


189 


29 


5 


do 


92 


. 66 


27 


5 




61.5 


48.5 


25 


2.5' 


do 


71 .5 


12.6 


28 


2.5- 


— — do — — 


63 


11.0 


28 


2. 5 


do 


59 


10.4 


- 28 


2. 5 


vertical 


58 , . . .... 


... 48 


28 




up 








. 2. 5 


same 


92.5 . 


66.5 


27 


2.5 


vertical 


61 


46 


26 




down 








2.5 


same 


;~ 58 ""' " 


63. 5 


26 



• *I"or a displacement of the flame front with velocity v 0 -b S , 
the quantity of substance consumed in the flame front per 
section is v obg p Q -rrrs grams. On the other hand at the 
flame front area S there burn each second v m S grams 
whence v obs Po nr 2 = v m S = v 0 p Q S. 
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The first method has the advantage that, besides per- 
mitting the computation of the flame-front .area S, it 
also permits the computation by the above ; formulas of the 
normal velocity of propagation v Q . The results of the 
computation are given in table I. 

We see that for tubes of different diameters for 
various shapes of the flame front and various observed 
velocities of the flame propagation the computation 
actually gives a constant (within the limits of obser- 
vation errors) normal velocity of propagation.* 

We have already shown that the uniform flame prop- 
agation in a tube is observed only at the start of the 
process over the first meter of the flame path. There- 
after the flame propagation, for reasons not yet under- 
stood, presents a nonuniform character. The flame, at 
times begins to move ahead with a very large velocity-; 
at times the velocity decreases and oven reverses its 
direction - that is, a kind of oscillation begins to 
take place. This is- clearly brought out, for example, 
in figure 23 (reference 9). It is also seen in figure 
3 where the regular flame front, after traveling 50 
centimeters along the tube, disappears suddenly and the 
picture becomes confused owing to the energetic flame 
oscillations. Sometimes the latter attains a very large 
amplitude as is seen in figure 5 (reference 9). In many 
cases there occurs another type of disturbance of the 
regularity of flame propagation in the tube as shown in 
figure 6, giving an example of the acceleration of the 
flame propagation and the occurrence of a detonation wave 
(taken on a rotating film, stoichiometric mixture of 
methane, photograph by Sokolik (reference 9)). As may be 
seen, the flame travels slowly at first with a constant 
velocity, then the velocity begins to increase gradually 
and finally at the distance S from the point of igni-. 
tion the flame' propagation enters a new. regime charac- 
terized by detonation. 



♦In wide horizontal tubes the convection in the direction 
of gravity disturbs the correct shape of the flame front, 
producing in it all kinds of nodes and breaks, a condition 
which renders impossible the measurement of the normal 
velocity of propagation. The lower photograph (fig. 12) 
of the article by Sokolik (reference 9) is an example 
of this type of flame propagation. 
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The detonating combustion (like the normal one) 
progresses at a constant rate with the difference, how- 
ever , that 'tfie letonating' velocity exceeds the normal 
comtmstion rate hundreds and thousands of times. Det-" 
onati.on does not arise in every mixture hut only in 
mixtures with a sufficiently large amount of heat liber- 
ation. The normal propagation is associated with flame 
propagation through heat conduct ivity, the flame front 
heating the neighboring cool layers of fresh gas to a 
high temperature, insuring a high reaction velocity? 
whereas in the detonating regime, the flame propagation 
occurs with the aid of the hydrodynamic shock wave. 
Over this wave front there is a sudden compression of the 
unburned gas to a very high pressure. As a result of 
adiabatic compression the temperature of the gas sharply 
rises and the mixture reacts suddenly. In our present 
paper we shall not consider these interesting phenomena 
which are important both from a theoretical and practical 
point of view. We shall only touch upon the problems of 
the fluctuating type of propagation which often arises 
to disturb the regularity of the normal propagation. 
There sometimes occur cases where the amplitude of these 
fluctuations is not large. In these cases, notwithstand- 
ing the presence of fluctuations, the average in time of 
the propagation velocity maintains the same value as the 
velocity of uniform propagation that occurs at the start 
of the process. According to Jost this may be considered 
as proof of the fact that for small fluctuations the flame 
front remains stable and maintains its shape. Relative to 
the unburned gas the flarne is propagated with the same 
constant velocity. The reason for the apparent irregular- 
ity of motion of the flame with reference to the observer 
appears to be the oscillation of the gas column in the 
tube as a whole. These oscillations originate from cer- 
tain unexplained hydrodynamic causes connected with the 
motions of the gases during the flame propagation. Small 
flame oscillations can be made out on figure 2 (successive 
flame fronts c ompres sed, then rarified). 

In those cases where the amplitude of the oscillation 
becomes large this oscillatory motion of the gas, by 
changing the distribution of the flow velocity over the 
tube cross section, leads in turn to a curvature of the 
flame front, increasing its area and thus also the mean 
rate of combustion. 

Flame acceleration preceding detonation, as Shelkin 
has shown, also appears as a result of the progressive 



8 



EA.CA Technical Memorandum No. 1026 



increase in area of the flame front and is not associated 
with any increase in the normal velocity. He showed that 
this acceleration is observed in the case where the gas 
motions, arising from the combustion, have a velocity 
exceeding the critical Reynolds velocity, that is, where 
turbulence may arise leading to an increase in the area 
of the combustion front and hence to an increase in the 
observed velocity of propagation. If the latter is so 
small that the critical velocity of the gas flow is not 
reached no detonation arises. If the normal velocity 
is sufficiently large nons imultaneous turbulent motions 
may arise in a smooth tube and gradually increase. Tor 
this reason observe at first constant normal velocity 
of propagation and only after the flame has traversed a . 
more or less long portion of the tube does nccel eration 
of the flame with detonation arise. The correctness of 
these views was demonstrated directly by Shelkin by 
observing the flame propagation in a tube with roughened 
walls (to obtain which he placed a wire spiral inside the 
tube). The turbulent motions were then expected to arise 
considerably earlier and this was actually found to be 
the case, the detonation being found to arise in the 
immediate neighborhood of the ignition sovirce at a dis- 
tance of a few centimeters of flame travel; that is, the 
distance from the source of ignition to the detonating 
point was shortened more than ten times as compared with 
a smooth tube. 

We thus see that every kind of disturbance in the 
observed velocity of propagation of the flame is connected 
either with fluctuation of the gas column or with an in- 
crease in the area of the combustion front due to hydro- 
dynamic causes, (for laminar flows by a change in the 
shape of the flame front and for turbulent flows by an 
overlaying of the flame front surface by sort of small 
ripples). The normal velocity of propagation remains, 
however, constant for all these regimes and detonation . 
appears princ ipally by the new type of the flame prop- 
agat ion. 

There are, however, more direct methods of deter- 
mining the velocity of flame propagation under conditions 
where the hydrodynamic state of motion of the gas is 
considerably simpler than for flame propagation along a 
tube. In the case of flame propagation in a spherical 
vessel with ignition at the center the entire arrangement 
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is symmetrical so that only radial flows are possible 
and there is n_o_. reason f or the spherical curvature of . 
the flame frqnt. 

Eadial flow arises from the condition that the gas 
burning at the flame front (thus enclosed within the 
flame sphere) has a combustion temperature many times 
exceeding the temperature' of the unburned gas and hence 
occupies a large volume. This forms a radial flow of 
unburned gas ahead of th'e flame front, the velocity of 
which, compounded with the normal propagation velocity 
(relative to the unburned gas) leads to an observed 
velocity T 0Ds considerably in excess of v Q . In the 

case where, for the spherical vessel a soap bubble is 
employed which expands with expansion of the gas, the 
process of combustion can be carried out at constant 
pressure and the computation of the observed velocity 
v obs as a Motion of the normal velocity v Q con- 
siderably simplified. 

Assuming, provisionally, an element of the flame 
front of unit area as stationary, we know that v m = v 0 p 0 
grams of unburned gas per second flows up to it, hence 
the same amount of burned gas is removed. If the density 
of the burned gas at the temperature of combustion is p, 
the velocity of the burned gas v multiplied by p 
should, according to the .law of conservation of matter, 
be equal to v m = v 0 p 0 , whence v Q p 0 = vp or v = v 0 p 0 /p, 

v o = v P/Po* But the velocity v with which the burned 

gas, moves relative to the flame front is equal to the 

velocity of the flame front relative to the burned gas, 

and since the burned gas under the test conditions is 

stationary relative to the observer, v = v . is the 

obs . 

observed flame propagation velocity. It may be measured 
by the two above-mentioned photographic methods, namely: 

(1) by a series of successive instantaneous pictures; 

(2) by photographing the flame motion (expansion of 
the flame sphere) on a rotating film. (See fig';. 16 of the 
above.-mentioned article by.So.ko.lik (reference 9) in which 
the soap buble at the beginning and end of the process is. 
indicated schematically). 

Thus ,• knowing ; v 0 ^, s ., we obtain the normal velocity 
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v obs P V r Q 3 

where T Q is the initial temperature of the gas and T x 

the combustion temperature, that is, the temperature of 
the burned gases. Since according to the Gay-Lussac law 
the volume v = 4/3- nr 3 (at constant pressure) is directly- 
proportional to the absolute temperature T the ratio 
Tq/Tj will be equal to the ratio r 0 3 /r 1 3 where r Q is 

the initial radius of the bubble and r x the radius of 
the biibble when all the gas in it is burned. (This radius 
can readily.be found from the photograph.) 

We thus obtain an elegant and very clear method from 
the theoretical point of view for the measurement of the 
normal velocity of propagation. If we operate with closed 
tubes or solid closed spherical . vessels then as a result 
of the expansion of the gases during combustion, the 
pressure as the flame travels, will increase. Since the 
propagation velocity may depend .on the pressure and since 
the hydrodynamic computation is, in this' case, complicated, 
the method of measuring the normal velocity under these 
conditions is very difficult. By employing long tubes 
or large spherical vessels, however, and restricting the 
measurement to the initial phase of the flame propagation 
(when the rise in pressure is not large) we may success- 
fully apply to the analysis of the results all the con- 
siderations and formulas which were obtained by us for 
tubes open at one end and for a soap bubble. 

The simplest method of all for measur iag the normal 
velocity is by means of the Bunsen burner, that is, the 
method of the stationary flame. Its disadvantage appears 
to be only the large expenditure of combustible gas. If 
the flame front in the tube were plane and hence propa- 
gated with the normal velocity v q we might, by blowing 
gas with velocity v 0 toward the flame, bring the flame 
to rest and thus obtain a stationary flame in the gas 
stream. The position of the flame, however, would not 
be stable because small fluctuations in the velocity of 
the stream would throw the flame forward or backward 
depending on whether the velocity of the stream was smaller 
or greater than v Q . If, in the axis of the tube, we 
place a point ignition source in the form of a hot point 
or spark it is possible to increase the velocity of the 
stream considerably above v Q and maintain a stable flame. 
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The flame front will then "be in the shape of a cone 
spreading in the direction of the flow, its vertex 
situated at the point of ignition and its base on the 
walls of the tube, The greater the velocity of the 
stream v the more extended i-s the cone. The cone 
angle, the angle between the axis (or the velocity of 
the stream) and the slant side, is such that v 0 = v sin cp 
It is readily seen that the quantity of the gas supplied 
"by the stream to a unit area of the combustion cone will 
be equal to p 0 v sin cp = p 0 v Q and this magnitude should 

as we know, be equal to the mass rate of combustion v m = 
v o Po* Prom this it is clear that the velocity of the 
stream may be considerably greater than v 0 without caus 
ing a break in the flame. If the ignition source is taken 
to be not a point on the axis but a heated wire ring place 
against the walls of the tube it will then be found with- 
out difficulty that the shape of the front should again be 
a cone having its base on the ignition ring. Both cases 
are sketched on figure 7, 

Tests, such as those described above, are much more 
easily and simply carried out with a Bunsen burner. At 
the rim of the burner stagnation points of the flow are 
formed where the velocity of the flow is small and there- 
fore when the burner is lighted these points play the 
part of the constant ignition source. The gas strea'm 
escaping from the turner (the burner usually operated 
with laminar flow) moves for a considerable distance in 
the form of a cylindrical column with the sane velocity 
v -that prevails within the burner tube. The conditions 
are thus obtained for the formation of a stationary 
conical flame front. Knowing the amount of fuel gas 
supplied and the cross-sectional area of the burner we 
determine the velocity y and from the cone angle cp 
find the normal velocity of propagation v Q - v sin cp. 
It is simpler to measure not the angle . cp but the height 
h of the cone. The angle cp , is then obtained from the 
equation tan cp = r/h (fig'. 8).'. 

The simple conical shape is obtained in the case 
where the velocity of the gas flow is uniform over the 
entire tube cross section. Actually with laminar flow 
of the gas along the tube the velocity distribution is 
parabolic, being a maximum at the center and falling off 
toward the walls of the tube. It will readily be seen 
that under these conditions the stationary flame front 
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will differ somewhat from the simple cone and will assume 
the cap shape shown in figure -9. This is, in fact, the 
shape assumed by the inner cone of a Bunsen burner. For 
this reason the measurement of the normal velocity, accord- 
ing to the previously given formulas from the height of the 
cone, involves a certain error. The more accurate method 
is therefore to photograph the cone and determine its area 
S from the photographs. Knowing the combustion area and 
the volume of gas per second V supplied to the burner 
the normal velocity is easily found from the formula 
v Q = V/S. 

In determining the propagation velocity by the 
Bunsen cone method it is often more convenient to make 
use of the arrangement with separated cones. 

The greatest number of measurements of the normal 
velocity refers to air and oxygen mixtures of carbon 
monoxide, hydrogen, methane or technical combustible 
gases of the type of illuminating gas. Less study has 
been devoted to the flames of the higher hydrocarbons. 
The experimental data thus almost exclusively refers 
to oxidizing reactions with free oxygen which for 
practical purposes are of most interest. The processes 
with other two-component systems (for example, combustion 
of hydrogen with chlorine and organic substances) and 
also with one-component systems, that is, the cases of 
flame propagation with exothermic decomposition of the 
gaseous substances, have received almost no study (except 
for the case of flame propagation in ozone). 

For two-component systems it is of great interest 
to establish the dependence of the normal velocity of 
propagation on the percent of combustible gas (CO, H a , 
CH4) for various ratios of the oxygen to the inert gas. 
On figures 10, 11, 12, 13 are. shown the corresponding 
curves for the mixtures CO, H 3 and CH 4 . The propa- 
gation velocity for CO and H 2 reaches a maximum, not 
for the stoichiometric ratios, of t he component s but for 
a marked excess of the combustible. Dilution of the 
mixture with inert gas lowers the velocity of propa- 
gation of the flame, (see figs. 11 and' 12) , a fact 
which is to be expected because the maximum temperature 
of combus tion is... 1Uie--r-eby reduced. The greater the heat 
capacity of the inert gas the more the combustion temper- 
ature is lowered and the more the flame propagation veloc- 
ity reduced. Tb'&s , if the mixture of methane and air is 
first diluted with carbon dioxide and then v/ith argon the 
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flame propagation velocity in the first case is about 
half that in. the second case. 

Preliminary heating of the mixture increases the 
velocity of flame propagation. On figure 14 are shown 
the results of Passauer for air mixtures of carbon 
monoxide containing 3.3 percent water vapor. The flame 
propagation in mixtures of carbon monoxide depends to 
a surprising extent on the water vapor content. The 
propagation velocity increases with increase in the 
percent moisture. This is shown in figure 15 according 
to the data of piock and King. A n addition of hydrogen 
and other substances containing hydrogen has the same 
effect. A mixture of carbon monoxide with air, or even 
oxygen that is perfectly dry and without hydrogen-con- 
taining substances, is altogether unsuitable for the 
flame propagation, that is, is not a combustion mixture. 
If certain substances, having a tendency to take up 
moisture, as for example, 0C14 , are added to a mixture 
of moist CO they strongly reduce the velocity of the 
flame propagation. These substances produce, however, 
the same action on the propagat ion velocity of hydrogen 
flames where the addition of moisture naturally has no 
effect . 

The problem of the variation of the propagation 
velocity with pressure has experimentally not been clari-. 
fied, the results of different authors being in sharp 
contradiction. Moreover, in tests with CO mixtures 
where the propagation velocity strongly depends on the 
amount of moisture, different results are naturally ob- 
tained depending on whether the tests are conducted at 
constant partial pressure of the vapor or at constant 
percent moisture. The latter case is of theoretical 
interest but in practice the first one is generally used. 
On analyzing the data of Ubbelohde, Zeldovich came to the 
conclusion that at constant percent moisture the normal 
propagation velocity of CO air mixtures does not depend 
on the pressure of the mixture. In CO oxygen mixtures 
the propagation velocity increases somewhat with the 
pressure, approximately, according to the law v Q « **/p. 

In mixtures of hydrogen with air, benzol and benzene the 
the propagation velocity apparently does not depend on 
the pressure. In methane mixtures the velocity drops 

with the pressure approximately as 1/ ./p. AH that can 
be affirmed with this state of experimentation is that 
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the velocity of the normal flame propagation in oxidizing 
reactions varies with the pressure according to the law 

v 0 « p n where n lies between 0 and -1/2. The corre- 
spondingmass rate of combustion v m = v Q p increases with 
the pressure according to the law p m where m lies 

between 1 and + l/2. Only in oxygen mixtures where the 
temperature of combustion is very high is there observed 
a slight increase in the propagation velocity v 0 with 

the pressure and a corresponding increase in the mass 
rate of combustion " v ffi with the pressure to a power some- 
what higher than one. 



2. ANALYSIS OF THE OLD THEORETICAL VIEWS 
OH FLAME PROPAGATION 



We have seen that the phenomena of slow flame pro- 
pagation all lead to the property of the flame front; 
namely, of a definite quantity v ffi grams of mixture 

burned per unit area of surface or, what amounts to the 
same thing, of being propagated relative to the nonburhing 
gas in. a direction normal to its surface with a propagation 
velocity v Q cm/sec. The observed motion of the flame is 

the result of the superposition of all kinds of hydro- 
dynamic motions upon this fundamental motion of the flame 
front. Thus, the theory of slow flame-front propagation 
is the theory of processes occurring in a thin combus- 
tion layer in the flame front. Since the curvature of 
the surface of the flame front is always, small in com- ■ 
parison with the thickness of the front (that is, the 
zone of reaction) we may with sufficient accuracy con- 
sider that all magnitudes (for example, temperature) 
within this zone are functions of a single coordinate, 
the direction of which is perpendicular to the element 
of front area (one-dimensional problem). Moreover, by 
associating the origin of coordinates with the flame 
front, that is., by studying the process in a system of 
coordinates moving in space with a given element of the 
flame front we may assume, that the distribution within 
the combustion and heating zone of all magnitudes of 
interest (for example, the temperature) does not change 
with time. The proplem thus resolves itself to the solu- 
tion of the one-dimensional stationary case. 
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For the computations bel ow, the following notation, 
will be used: 

T Q , Ti respectively, the initial temperature of 

the unturned gas, the variable temperature 
in the process of combustion at various 
points of the flame zone, and the maximum 
combustion temperature without taking heat 
losses into account. All temperatures are 
.expressed in degrees on the absolute scale- 
's the difference T x - T 

a Q , a the number of molecules per unit volume of 

the combustible in the initial cold lixc 
ture and at various points of the combustion 
zone, respectively, in the case where a 
decomposition reaction occurs (C1 S 0, 0 3 , 
explosive substance) 

n, n Q total number of molecules in the initial mix- 

ture at temperatures T and T Q >-- respectively 

A coefficient of heat conductivity of the gas. A 

increases with the temperature 

A" coefficient of the heat conductivity of the 

products of combustion at temperature T x 

Cp specific heat of the gas at constant pressure. 

Cp increases with the temperature. At first 

we shall assume Cp independent of the 

temperature and equal to its average value 
between To and Ti 

c£ c specific heat of the products of combustion at 

the temperature" Ti 

.H, H 0 i H x heat contents of unit mass of. the gas corre- 
sponding to temperatures T, T Q , Ti 

p Q » p density of the gas at temperatures T 0 and T, 

respectively 

p* density of the products of combustion at 

temperature Tx 
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normal velocity of flame propagation in centi- 
meters per second, that is, the velocity of 
motion of the flame front relative to the ua- ■ 
burned gas along the normal to the front sur- 
face at a given place (or, what amounts to 
the same thing, the component along the normal 
to the front of the velocity of motion of the 
cold gas relative to the flame front) 

the velocity of the gas relative to the flame 
within the combustion zone at temperature T 

the mass rate of combustion, that is, the number 
of grams of the gas mixture burned per 1 cm 3 of 
the flame front 

the coefficients of diffusion of the gases at 
■ temperature T and T Qj respectively 

the fraction of the weight of the combustible in 
the initial mixture or equival en tly the number 
of grams of combustible per gram of the mixture. 
In the case of two-component systems, for example, 
for oxidizing reactions M 0 will be the fraction 
of the fuel that oxidizes in the process of com- 
bustion. M Q is always < 1 

the fraction of combustible per gram of the mix- 
ture at various points of the combustion zone 
M .< 1 

the Avogadro number = 6xl0 23 molecules 

the heat of reaction of a single molecule of the 
combustible and in the case of two-component 
systems, for example, the oxidation of CO to 
C0 3 , the heat of oxidation of CO referred to 
one molectile of GO 

the heat of combustion of 1 gram-molecule of com- 
bustible (the heat of decomposition of i gram- 
molecule C1 2 0, the heat of combustion of 1 
gram-molecule of CO etc.) 

the heat of combustion of 1 gram of the initial 

mixture (for the combustion of a single or two- 
component system diluted with an inert gas; L 
refers to 1 gram of the entire mixture including 
the inert gas 
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the speed of reaction expressed "by the number of 
reacting molecules of the combustible per unit 
volume per second. Jor the reactions of zero, 

first and second order w - Se - ^/^, w = 'ka.e~ £j ^'~ 



and w = ka e 



-B/ET 



, respectively 
S energy of activation of the reaction 
E gas constant 

(j, molecular weight of the combustible 

Z impact coefficient, that is, the number of. impacts 
per unit time per unit volume between two mole- 
cules if the unit volume contained just these 
two molecules 

I free path between impacts 

u velocity of the heat motion of the molecules 

<y effective diameter of the molecules on impact 

m mass of a molecule of the mixture 

m mean mass of a molecule of the mixture 

[The following simple relations between the above 
magnitudes will be required for what follows: 



M oPo = 



a 



Po 
- toil 



?o- 



M 0 ft M 0 q»M M 0 q' 



m 



= H, - H, 



0, = s<i«- 



±0 



Por constant specific heat L = Cp (T-y - T 0 ) 
For variable specific heat H =' / c p dl ; E 0 ~ / c T d T ; 

an 



> (a) 



0 



Hi = / c p dT; c p = ~~y D = 1 lu; A = ^ *upc p ; 



cr 



1 _ 

cjsun Zn 



am „ a Q m 
= M; — — 

P Po 



M 0 1 3STm = n ; nm = p 
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Finally we present various possible definitions of 
the fuel concentration and the temperature in the com- 
bustion zone. 

1. a the number pf molecules of the combustible 
per unit volume at different points of the combustion 
zone. This is the concentration we shall mainly use. 

2. — = C relative molecular concentration, that 
n 

is, the fraction of the molecules of the co mbus t ibl e in 
the mixture at various points of the combustion zone. 
■jj x 100 percent relative concentration. 

3. M the weight concentration of the combustible 
during combustion (if the molecular weights of all the 
gases composing the mixture and of their products of 
combustion are approximately equal to each other 

M = a/n = C) , B = M/M 0 the relative weight concentra- 
tion of the combustible. If the molecular weights of 
all the gases in mixture are the same then 

M a a 0 C 
,M D n " n 0 Co 

4. 3? the absolute temperature at various points 
of the combustion zone. 

5. 3? - T Q rise in temperature a.bove the initial. 

6. 57 = (3? - 3? 0 )/ ( 37 x - 37 0 ) relative rise in tem- 
perature above the initial. 

In the flame front within the heating and reaction 
zone the temperature rises from the initial temperature 
T Q of the cold gas to the temperature 37 x , the maximum 

combustion temperature measured in thousands of degrees. 

Considering the temperature distribution within the 
combustion and heating zones it is not difficult to arrive 
at the following differential equation 



(1) 
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In the -.case where X / c may he considered a constant 

magnitude within the temperature range T Q to T a (that 
is, in an interval of the order of 1500° - 2500°) the 
equation may be written in the form 



iL All _ dH + w , = 0 (lf) 



c p dx* 5 111 dx 

Assuming c constant within the range T Q to T x we 
arrive at the equation 



0 

+ — 

dx s dx 



d s T dT wQtW ^ 

- a — + = 0 (2) 



where 



a - x 



then we shall derive this equation. When in the gaseous 
phase heat is evolved and the hea.t is transferred from 
some points of space to others then (even in the absence 
of forced motion and convection under the action of 
gravitational forces)) there must occur a process of 
transfer of individual elements of the gas as a result 
of the heat expansion. Thus even in this case, and 
all the more in the general case, to the heat conduc- 
tivity must be added heat transfer arising from the 
gas motion. We shall denote' the linear velocity of 
the latter by V (a vector which varies in direction 
and magnitude at different points of the space occupied 
by the gas). In this case the flow of heat per unit 
time through an area of 1 cm 3 perpendicular to the 
direction of 7 will be pHT where pH is the heat 
content of a unit volume of the gas. Moreover heat 
will be transferred by the heat conductivity, the amount 
being, as is known, X grad T. The total heat will be 
equal to 

Q, = -> X grad T + pHV 

If the reaction proceeds with velocity w then 
each second wQ« cal of heat energy are given out per 



I 



20 NACA Technical Memorandum Ho. 1026 

unit volume. Hence the change in heat content per 
unit volume per unit time Q (pH)/dt will be equal to 

a ( P H) = _ div ^ + _ div \ grad rp _ div ( p HV) + wQ« 

dt 

or making use of the equation of the conservation of 
matter 

- |£ = div (pV) 

and of the mathematical relations 

SiEal = p-P+ H§£ and div (pHV) = H div ( p V) + pV grad H 
d t d t o t 

we may rewrite the equation of heat distribution in the 
form 

p |^ = div X grad T - pV grad H+ w!j» 
Using the relations c p = dE/dT and c grad T = grad H 



♦Strictly speaking this 'is not so becattse in the case of 
chemical reaction where the composition of the mixture varies 
in time and space the heat content H varies not only with • 
the temperature but with change in M. If Cj and c 2 are 
the specific heats at constant pressure of the initial sub- 
stances and final products , "respect iveiy , then 

h-h 0 = /[.,»♦ (!-«..] « or || =(|f) M §| + (§f\ ||; 

To T 

(I) M = c p - + c 1 -^ (5), - /Wo> d * 

' To 

From what follows we shall see that under the flame Con- 
di t i ons 

d x M 0 dx 

where Ti is the maximum temperature developed in com- 
bustion, whence (Footnote continued on page 21) 
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we may write 

div X grad T = div -~- grad H 

°P 

and rewrite the fundamental equation in the form 

3 H X 
p — - <iiv — grad H - p Y grad H + wQ' 
3t c p 

- — -' — ■ ■ 

* (Footnote continued from page 20) 

/3H\ 3M . /3H\ ST _ 
MM^i Qx ' VaT-i ox 

T 

M 0 ^ (c, - c 3 )dT 



_ To Mp (c a ~ c» ) ~T n ) 

(T a - T 0 ) c p = T, - T 0 c p 

where c x and c 3 are the mean values of the specific 
heats in the temperature range T 0 to T. Chemical trans- 
formation, as we shall see 'below, only proceeds at high 
temperatures for which the difference c 3 - c t is always 
small. It is not difficult to prove from tabulated data 
that even at temperatures of some hundreds of degrees the 
difference c 3 - c x for the majority of cases of interest 
to us is less than 10 percent of' the value of c^. The 

T — IP " 
magnitude -°- M 0 is always less than unity. We may, 

Tj - T 0 

therefore, neglect the term -?2! in the expression for 

3M ox 

|| and consider that |§ = (||) m -g = Cp || or, general* 

ized, c p grad T = grad H as was to he proved. Tor the 

same reason the heat Q' or Q given off by the chemical 
transformation is ,con stant ,^ not depending, on_ the .tempera- 
ture, at which the transformation occurs. 

We have, in fact. 

^T = ^T Q + ^ / (°i ~ c s) dt ~ $T 0 

in the case where Qrp. is sufficiently large as is 
always true in combustion reactions. 
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or in the one-dimensional case where all the magnitudes 
depend on a single coordinate x and the velocity vector 
forms with the x-axis an angle cp (for example, V gr.ad 
H = V cos cp dH/ dx) 

p = — I 2 - PV coscp j£ + WQ« ( 2 .) 
K dt dx Cp ox r r ox ^ > ' 

In applying to the flame zone where as we have seen the 
problem reduces to the one-dimensional and stationary case 
(dE/dt = 0) and where pV cos cp is no other than the 
weight of the burning gas mixture per second per unit area 
of the flame front, that is, pV cos cp = , we obtain' the 
equation given above. 

Previous to the work of lewis and von Elbe (reference 
2) and Zeldovieh and Jv ank-Kamenet sky , partly due to 
insufficient understanding of the process, partly for the 
sake of mathematical simplification of the problem, every 
author introduced into the theory a certain constant, 
namely, the self-ignition temperature of the given mixture. 
It was supposed that the fresh gas, on approaching the 
combustion zone was heated up by the heat given up by the 
latter to the self-ignition temperature T± after which 
the inflammation took place, that is, an intense reaction 
at a constant rate w was initiated in the mixture until 
all the substance was consumed and the maximum tempera- 
ture attained. As we have already seen above(see sec. 
1) the self-ignition temperature is not a constant but 
depends on the test, conditions, chiefly on the conditions 
under which heat is liberated and on the induction period 
of the zone of inflammation. Moreover, the- a s sumpt ion 
of a constant speed of reaction, on increasing the temper- 
ature from Hi to I a likewise does not bear criticism. 
Jot this reason the introduction of the concept of self- 
ignition temperature in the problem of the velocity of 
flame propagation is very reminiscent of the views with 
regard to the "reaction point" in the s elf- ignit i on theory 
prevalent before van 1 1 Hoff. 

Wo shall try, holding to the premises of the previous 
authors, to present these old theories in their strictest 
and most logical form without following, in too much de- 
tail, the often confused and not always correct concepts 
which they employed. We stop, to consider these' old theories 
based on doubtful premises not only because, until now, 
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they have been used by scientists and technicians, hut 
also because their clarification will be of value in 
presenting, in the next section, the new more rigorous 
thermal theory of flame propagation, the pr inc iples . of 
which are explained in the paper by Zeldovich and Frank- 
Kamenetsky (reference 3). lor the solution of the prob- 
lem we shall make use of equation (2). 

Starting from the above naive assumptions we are 
to formulate the boundary conditions under which equation 
(2) is to be integrated. We shall lay off the temperature 
of the gas on the coordinate axis and on the axis of ab- 
scissas the distance x from that point in the zone at 
which the inflammation temperature T^ is attained (fig. 
16). The x-axis is directed from the cold initial gas 
to the products of combustion. In our assumed system of 
coordinates moving with the flame front the latter is 
stationary and the fresh gas moves along the x-axis with 
the velocity v Q toward the flame front. At x = -co 
the cold gas has the initial temperature T Q . Then in 
the heating zone and moving in the direction of positive 
x the gas is heated by conduction up to the temperature 

(region I). The gas then starts to react and continues 
to be heated up (region II). After traversing the distance 
d from the origin of coordinates the weight of the initial 
gas has completed its reaction and the temperature attains 
the maximum combustion temperature T a which in the absence 
of heat losses remains constant further along (region III). 

As we have already shown, the reaction velocity w 
at the self-ignition temperature and above has tacitly or 
explicitly been assumed constant by the previous authors. 
It is immaterial whether to assume as constant the number 
of transformed molecules w per unit volume per unit time 
independent of the gas density which varies as the gas is 
heated or the time taken to transform each given molecule, 
that is, to consider w/p = constant or w proportional 
to p (which changes with the heating). The second 
assumption leads in the final expression for the propa- 
gation velo city t o a value less than the first in the 
ratio 1: J T Q / T ^ , that is, approximately 




The first assumption is equivalent to a reaction of zero 
order, the second to one of first order with respect to the 
density or pressure. As we shall see below, for a reaction 
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of. the first order it is necessary to take into account 
the diffusion in the zone as was not done in the older 
theories. Only for this reason, in fact, is. the assumption 
of zero order of the reaction in the older theories de- 
void of inner contradictions and we shall, therefore , 
stop to consider it. 

Let us consider the motion of a small element of 
volume dcu 0 of the initial cold gas at temperature T Q 
containing a 0 dtu Q molecules of combust ible capable of 

chemical transformation. This element of volume gradually 
heats up as it moves from x = - » in the direction of 
positive x the volume of the element increasing due to 
heat expansion so that at a certain distance x and at 
a corresponding temperature T the volume of the element 
will be 

m 

du) = — duj 
T 0 0 

For T > Tj_ a chemical reaction in the element of volume 
do> will take place with velocity wduu, where w by 
assumption is a. constant magnitude. In t-he interval dx 
between x and x + dx the element of volume remains 
a certain time dt evidently equal to dx/v where y 
is the velocity of motion of the volume element at the 
given point. If the cold gas moves with velocity v 0 
then on increasing the temperature the gas due to e x- 
pansion will move more rapidly, namely, with velocity 




so that dt = dx/v 0 x T Q /.T 

The number of reacting molecules in the element 
dx is evidently equal to 

dx Tn w 
w dU) dt = w dO) = dx cto 

v o *i v o 0 

that is, in each element dx the same number of mole- 
cules will react. Therefore, the total number of reacting 
molecules in the distance between x = 0 and x = d of 
the volume element, will be wd/v Q x du) Q . But by defini- 
tion all molecules a du) Q of the combustible in the 
element of volume at She distance d are reacting, whence 
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^ = a 0 or d = ^Zo (3) 
o w 

We thus obtain the required value of the distance d. 

In integrating equation (2) vie divide the integrating 
range into two subint ervals . In the first one where the 
temperature rises from T = T 0 to T = T± and x varies 
from - oo to 0 no reaction, according to the assumptions 
made, occurs and hence w = 0. In the second interval 
in which T increases from T = T^ to 3? = T x and x 
from 0 to x = d the speed of the reaction, by the 
assumption made, is constant. In the third region from 
x = d to +oo and w = 0 the temperature is constant 
and equal to Tj and therefore the solution is obvious. 
Thus for region I 



^1 - a A* = 0 (4) 



dx 3 dx 
where 

for x = - co, T = T 0 and for x = 0 , T = T i (5) 
In region II 

d 2 T dT wO» 

dx 3 dx A 1 ' 

where 

for x = 0 T = 1^ and for x = d T = Tj (5») 

As may be seen these conditions are quite sufficient for 
the solution of the problem. Besides these, however, 
there is another condition, namely, the condition of 

continuity of the heat flow q = - X — where the two 

dx . 

regions meet, that is, the condition 



MxA=0,I = Wx=0, 



II 
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This added condition determines the value of the param- 
eter. v Q or v , that is, the normal velocity of the 
flame propagation. W e find the equation for the deter- 
mination of v m or v Q = v m / p to "be 



T ± - T 0 1 - e~t 



T 1 - T 0 



(6) 



where 



e _ , T m c p d . a o v o S Po c p - v o 3 po c pT 

£ = ad = A "IT A " A (6t) 
and the magnitude a 0 /w is denoted by T. 

We shall show how the solution is found. . Accord- 
ing to formulas (A) and (3) 

wQ' wL(a (Tj - T Q ) T x - T 



o 

= a 



A M Q Ito dA d 

Substittxting the above in equation (4) and integrating 
under conditions (5) vie find for region I (x from- <x> to 0) 



(T - T 0 ) j = (Ti - T Q )e ax (?) 

and for region II (x from 0 to d) 

T, - T 0 ad ax 

( T - To)lI = _ x + ( Ti - T 0 ) 6 ^ - (7 .) 

forming (dT/dx)j and (dT/dx)jj for x = 0 and equating 
them, we obtain 

10 d e ad - 1 

whence v/e obtain formula (6). 
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When T i is near T 0 , that is, the self-ignition 
temperature is low 



« 1 

*i " *o 



This will he the case when e - ^ is small or | is 
large. Then 



Ti - T 0 1 



*i ~ T o " ^ ~ ad ~ v 0 2p 0 c p f 



that is , 




*i - *o 



p T Tj - T 0 (8) 



P^o 



When T i is near 1 1 , that is, the self-ignition 

temperature is high, then ( T i -T Q ) / ( T X -T ) is approxi- 

mately 1. This will be the case when g is small and 

e- I in formula (6) may be developed into a series. 
Taking the first three terms of the series, we obtain 



.2 



*l - T 0 | 2 ~ 3 



(we neglect the term g$/6 £ a2d 3 /6). Thus 
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The formtila 




c pPo T Ti - T Q 



(10) 



which is arrived at by most of the previous authors is 
an approximation, as is readily shown, of the expression 
of the general formula (6) and, generally speaking, not 
accurate. 

In what follows we shall "be especially interested 
in the case where Tj_ is near the value T. . In this 
case formula (9) is applicable. We shall, estimate the 
range of application of this formula. We set £=ad=0.5. 

£ s 0 25 

Then neglecting the term = — ^ — = 0.04 we make in 

b b 

comparison with |/2 = 0.25 an error of 16 percent in 
the determination of the magnitude T x - t//T^ - T Q or 

8 percent in the determination of the velocity v . 

Thus the formula is applicable with the required 

accuracy up to ^ 0.25. If Ta - T n = 2000 

T- - T 
x i x o 

the formula is applicable up to T x - T Q = 1500° or 
Tj. - T x = 500°. 

In what follows we shall find of importance another 
approximate method of solution for the case when T± is 
near T x . This method leads to the same result as formula 
(9) and hence gives the same degree of accuracy. It 
differs, however, in this respect, that it may with equal 
accuracy be applied to cases where the reaction velocity 
w is not constant but varies, for example, with the 
temperature. This method was proposed b;r Zeldovich and 
Prank-Kamenet sky and is based on the assumption that for 
Tj_ near in value to Ti in equation (4') for region II, 
it is possible to neglect the second term and v/rite the 
system in the form: 



dx 3 



dT 

+ a — 

dx 



0 



(region I) 



(11) 



dg; 

dx 2 



+ 




0 



(region II) 



(11') 
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The validity of this will he clear from a consideration 
of the complete equation (4 ' ) if there is substituted 

in it for the magnitude wQ'/X its equal a ~„ ,. ? 0 . 

a 

Integrating the equation with respect to x from 0 
to d we obtain 

mi 4- ^ dx 4) 



= -a (^-Ti) + a (T x - T 0 ) 

Since T x - T x « T x - T c the first' term on the right 
side of the equation may be neglected, which means that 
we may neglect the term a dT/dx in equation (4')t that 
is, consider the term a dT/dx small by comparison with 
wQ'/JS and dST/dx 3 .* 

The physical significance of neglecting the term 
is the following; In region I the convective heat flow 
expressed by the term cz(T x -T D ) is directed toward the 
positive x values and the flow due to conductivity, 
expressed by the term dT/dx, is oppositively directed, 
both flows compensating one another and thus establish- 
ing a stationary state. At the boundary of the region 
the flow of heat by conduction, expressed by tho term 
(dT/dx) Q , becomes equal to a(Ti~T 0 ) - a( T x -T 0 ), that is, 
practically all the heat in the reaction zone is used in 
heating the fresh gas. Since the heating in the reaction 
zone (from temperature Tj_-T 1 ) is very small the heat 
in region II is conducted away without contributing 
practically to the further heating of tho gas (although 
the maximum reaction temperature is there attained). 



* It may appear strange here why we may consider the 
magnitude a dT/dx small by comparison with d2Q}/&x 2 
in region II while in region I these magnitudes are 
obviously of the same order. This is explained by the 
fact that on the boundary of the region at point x= 0 
the second derivative undergoes a discontinuity and for 
T x - T x « T x - T Q in region II it is much greater than 
the second derivative in region I while the first deriv- 
ative in region I while the first derivatives in both 
regions are of the same order. 
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Hence in this region equation (11') is approximately 
satisfied. 



Integrating equation (11) there is obtained 



v oPo L 



since "by assumption is near Tj_ 

Integrating equation (ll 1 ) there is obtained 



- (m\ ♦ civ gd„ = /V / - « 

Ti 

Equating the two expressions obtained for (d?/dx) 0 , we 
have the equation for the determination of a or v Q : 



— ^=y^X" f w d ^ or v o = -~ y^Aq' y°w dT (12) 
Writing the above in somewhat different form 



v 



* /gAQ.'ao 1_ 



^ c p Po V 1 ! i O'' 



where 



T 



I B — / w dl (13) 

a o 



Since 



a0 Qf = _ = Po „ p = PoL = p 0 c (T -T 0 ) 

Iff u Po m 0 h Po 
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v o= / — — : - i^t < 14 > 

Thus expressed the formula is applicable also for w 

w T» — T i 

varying with x. If w = const,, then I = — (Ti-Ti)^ 



o 



and formula (14) leads to formula. (9) earlier obtained by 
another method. 'It is thus seen that the assumed approx- 
imate method of integration has the same accuracy as the 
expansion of the exact solution (6) in a series to the 
third term inclusive. 



3. NEW THEORY OF FLAME PROPAGATION 



In the old theories, as we have seen, two assumptions 
were made: 

1. the mixture starts to react on the attainment 
of the temperature Tj_ and 

2. the reaction proceeds at a constant speed w. 

We know that the reaction speed as a function of 
the temperature is expressed by the formula- 



where E is the activation energy, a magnitude charac- 
teristic of a given fuel. For substances suitable for 
combustion this magnitude is usually high and has the 
value of 25 to 80,000 calories. We should ,. therefore , 
expect a very rapid rate of increa.se in the speed of 
reaction w with the temperature and, with the coordi- 
nate x in the flame. Hence, the assumption of non- 
dependence of the reaction speed on the temperature in 
the t emper atur.e interval., JLi, . t,o. T t is . in., no ..way jus t if ied, 
Thus if E * 60,000 calories, T x = 2000° and T ± = 1500°, 

the speed of reaction increases 120 times in passing from 
Ti - T x . If E = 30,000 calories the speed of reaction 
for the same conditions increases almost 10 times. In 
the case where T x = 2000° and Tj = 1700° the speed in- 
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creases almost 15 times for B = 60,000 and 3 times for 
E = 25,000 calories. 

With such a sharp decrease in the reaction speed 
with decrease in the temperature there is no necessity 
at all for introducing the self-ignition temperature. 
As a matter of fact this temperature was needed only 
for establishing a limit below which no reaction occurs. 
In actual cases, however, when the reaction speed drops 
very sharply with lowering in the temperature according 
to the law w = Se E / aT the latter itself determines the 
range within vrhich practically no reaction occurs. For 
this reason we may introduce the notation T^ = T' but 
no longer as a physical magnitude but as a mathematical 
device for approximate computation, namely, as that tem- 
perature below which there is practically no reaction. 
It is necessary . for us to show that the final result, 
that is, the computed velocity of propagation will be 
practically independent of the value chosen for this 
fictitious magnitude. It is to Zeldovich primarily that 
we owe this restatement of the problem.* 

♦Some perplexity is felt by many to be occasioned by the 
following problem. It is well known that the speed of re- 
action of most fuels already reaches a very considerable 
magnitude at temperatures of the order of 1000°. Never- 
theless at a flame temperature of 2000° we neglect the 
reaction speed at 1000° and consider that the entire re- 
action proceeds only at temperatures near 2000°. The 
matter is here very simply explained. The velocity of 
flame propagation, as we have seen, is proportional to 
the square root of the reaction speed. Hence the greater 
the latter the greater the flame propagation velocity. 
Jor this reason the latter will be fundamentally deter- 
mined by the maximum possible reaction speed, that is, 
the speed near the temperature T x . The greater this 
speed the. smaller the time interval within which the gas 
remains in the reaction zone. This interval is so small 
that in those parts of the zone where the temperature is 
considerably below T x the reaction has no time to pro- 
ceed at all, although for large time intervals and in 
closed vessels the reaction speed is sufficiently great 
at the corresponding temperature. Thus the smallness 
of the reaction speed in the zone is determined as its 
smallness with respect to the reaction speed near the 
maximum combustion temperature. 
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, . Under the- new assumptions the fundamental equations 
of heat distribution (l) and (2) remain unchanged except 
that the speed of reaction w will... depend on the temper- 
ature according to the law w = Se and in the limit- 
ing cases instead of there will enter the above men- 
tioned temperature 27'. In the case where the energy of 
activation E is sufficiently large, namely, when 
E/HIj^ » 1 the temperature I' below which the speed of 
reaction may he neglected will lie near T x (which cor- 
responds to the case where (T x — Ti)/(T X — T 0 ) « l). 

Denoting by 9 the difference between T x and any 
temperature T within the reaction zone (between T x 

and T' ) wo,nay a ( since 9 « T x ) set o • E / Sx ' equal to 
-VSTi -E9/RT X 2 1 



e , e . Thus 



w x = Se e ( 15 ) 

where w x is a constant magnitude equal to the reaction 
speed at the maximum combustion temperature T x . We see 
that the reaction speed drops e times for 9 = RT X 2 /E, 
7.5 times for 9 = HT x 3 /E x- 2 , and so forth. lor E = 
60,000 calories and T x = 2000° P.T X 2 /E £ 125°, Hence 
for 9 = 250° or T * = 1750° the reaction speed w drops 
almost 8 times compared with its maximum value w x for 
1 = 3?! = 2000°. Thus, if we restrict ourselves to this 
degree of accuracy and neglect the reaction speed at 
T < 1750° we commit an error not exceeding 12 percent in 
the determination of v Q 2 and still less in the determina- 
tion of v 0 . E.or such value of . 9 (= 250°) we may: l) 

cons id « 1 and make use of the equation (15) and 

2) due to the smallnoss of T x — T f by comparison with 

T x — T Q in the solution of the problem make use of the 

approximate system of equat ions ( 11 ) assuming in them 

-E9/RT, 2 
w = w x e ' * 

We then obtain a solution in the form of equation. 
(14) where in the given case 

T E9 0 E9 Pi'" 



E 



E8 

rtT 2 " 




ET X 2 _ w 



dT -- 

a o 



r - etps Wl EIx s r e 
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where $ % = 18 X /E!E :L 8 

If we choose 6i .eq.ual to 2RTi s /E, that is, neglect 
the reaction speed when it is 1/8 as large as the maximum 
(wj_) then p = 2; if we assume 9 = 3RT 1 2 /E , that is, neg- 
lect the reaction speed w when it is 1/2 as large as 
w : then p = 3; if 6,^ = 4ET 1 3 /E, that is, w'/wj = 50; 
then p, = 4. 

j = / e~P dp = 1 - e~Pi 




If p.! = 2, j = 0.87; if 0 t = 3, J = 0.95; if p a = 4, 
j = 0.98. Thus we see that the choice of 63. provided 
61^ 2ST 1 3 /B has practically no effect on the magnitude 
j and hence on the speed of reaction. We may thus with 
an accuracy sufficient for our purposes assume in every 
case 3=1. Otherwise expressed we may everywhere in 

formulas (14) and (12) set / w d I equal to / w dT. 

T « 0 

This very important result permits us to make use of the 
auxiliary magnitude T ' without assigning to it any 
fundamental significance. In the final result it drops 
out. Thus 

_ ll aT i S . Se~ E / HTl ETi s . 
a 0 E J a 0 E 

whence according to formula (14) 



v 

o 



= / 2* Se-S/STi ETjg , 4 Se^lEM (l6 ) 
ycppo a 0 (T 1 -T 0 ) 1 J p 0 a 0 L E 



The above equation agrees with formula (9) if the 
reaction speed w is taken to be w x = Se~ E / RT i at the 
maximum combustion temperature and the difference T - Tj. 
the magnitude RT^/E. In this manner Zeldovich succeeded 
in ejecting from the formula for the propagation velocity 
the doubtful constant T i and the ill-defined magnitude 
T, the reaction time. 

Even in the above form, however, the theory may not 
yet be considered satisfactory because the magnitude S_ 
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in it was assumed by us to "be constant (reaction of zero 
order). Actually, for all homogeneous reactions S depends 
in some way on the concentration of the reacting substances. 
In the case where the combustion represents the reaction of 
a single substance (decomposition of azomethane, C1 3 0, 
etc.) the magnitude S is some function of the number of 
molecules of the combustible in a unit volume of the mixture 
at a given temperature. In the case of monomolecular re- 
actions S = ka, for bimolecular reactions S = ka 3 . When 
the initial mixture consists of two reacting kinds of mole- 
cules S may be a function of the concentration of both 
kinds. For bimolecular reactions in this case S = kab. 
The number of molecules a and b are connected with each 
other by the laws of stoichiometry and diffusion.* 

The number of molecules of combustible a in the 
reaction zone decreases as the reaction progresses. More- 
over the reaction products diffuse from the reaction zone 
into the unburned gas and this also decreases the value of 
a. The dependence of a on the coordinate x in the 
flame region may be found from the simultaneous consider- 
ation of the equations of diffusion and heat conduction 
(as was done by Zeldovich and Jrank-Kamenetsky) leading 
to similarity of the temperature and the concentration 
fields. The same similar ity was postulated considerably 
earlier by Lewis. In the clearest and simplest form this 
was done by Zeldovich and Frank-Kamenetsky and these authors 
will be followed in the presentation below. There arises here 
however, a very difficult problem with regard to writing down 
the equations of diffusion and propagation where a tempera- 
ture gradient exists. Zeldovich writes the expression for 
the diffusive flow in the usual form D da/dx where D is 
the coefficient of diffusion of the molecules a into the 
reaction products. If the special phenomena of thermodiffu- 
sion are not taken into account it is evident that for equal- 
ity of the partial pressure of the molecules .a- there can 
be no diffusion flow, whereas in the presence of a temper- 
ature gradient the number of molecules in' a unit volume 
will be different at different points.** In this case it 

♦We exclude, for the present, from consideration the case of 
t.w t o„ ^component s and take, it up,, below .in section 3 in the exampl 
of the oxidation of CO. 

**0bviously this must be understood as follows. In the field 
of variable temperature there exists, as usual, the flow D 
da/dx which, however, is compensated by the reverse mass 
flow va due to the tendency of the system to balance its 
pressure. (Footnote continued on p. 36) 
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is necessary to write the equation of the diffusive flow 
expressing it through partial pressure. 

When the number of molecules does not vary during 
the reaction (we at first limit ourselves to this case) 
we may, instead of the partial pressure, make use of the 
relation between the number of molecules and density of 
the gas and write the expression for the diffusive flow 
in the form Dp d(a/p)/dx. In this case the change in 
the number of molecules of the combustible with motion 
of the gas across the flame front (along the coordinate 
x) is given by the equation 



d 

dx 



j a 
d P 



vp 



_E 

dx 



- w = 0 



(17) 



the derivation of 
of equation ( 1) . . 



which is entirely analogous to that 



** (footnote continued from page 35) Thus if the number 
of molecules a* in the equilibrium state is distributed 



Vp = a o/po 



as the gas density, that is, a ! 
a 0 p/p 0 then the equality of the diffusion 
may be written in the form of the equation 



or 

and 



a* : 

mas s 



f 1 ows 



D~— - va* = 0 

dx . 



whence the velocity of the mass flow is 



D da* D dp 
a* dx o dx 



When the number of molecules a changes' in space 
not in correspondence with the formula a/p = a 0 /p 0 we 

have to deal with true diffusive flow which is composed 
of the absolute flow D da/dx and the oppositely directed 
mass flow which balances the pressure and is equal to 
va where v, as before, is equal to D/p x dp/dx. Thus 
the complete diffusive flow is equal to 
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Since Dp does not depend on the density of the gas 
and hence to a first approximation also on the temperature 
(with the same accuracy as K/c ) and since vp = v m we 

may vrrite the equation of diffusion of the molecules of 
combustible in the flame front in the following form: 



ill 

Dp —4- 

dx 



d a 
dx 



w = 0. 



(18) 



The first term expresses the mutual diffusion of the 
molecules of the combustible and the reaction products 
under the action of the difference of concentration ahead 
of and behind the combustion front. The second term ex- 
presses the flow of the molecules of the combustible 
under the action of the gas flow and the third their dis- 
appearance due to the reaction. 

For x =- 00 (a/p)_ co = a 0 /p 0 . We introduce a new variable 
a = a 0 /p 0 - a/p and in equation (2) a new variable 
Cp . . 

6 = qf vT-T 0 ), The two equations then assume the following 
form : 



Cp dx 3 



- v 



m 



d_6 

dx 



+ w = 0 



(19) 



Dp 



.3 

d q _ 
dx 3 



'm 



da 
dx 



+ w = 0 



(19') 



Prom the kinetic theory of gases we know that to a 
first approximation Dp = A/ C p» The two equations then 
become identical. x 

As we already know the boundary conditions for equation 
(19) will be 



( 9> = ©T (9)^ = 2* (T-.-T ) = — = M 

Jor equation (18) we have the boundary conditions 
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(?) 



= ±9., i. e. (a)_oo = 0 and (a,) +0o = — ft 



P/-a> Po 



But according to formula (A) 



apu, 1 or ao MpH 
lra oPo Po »* 



that is , 



(e)-co = (c0_co = 0; (9) +CB = (a). 



Po 



If the equation and the "boundary conditions for 6 
and a coincide, that is, a = 6 in the entire interval 



a c p T 



_ CpT 0 



or 



c p T + ^ = c p T 0 



iQ, 1 



(20) 



H = CpT is the store of thermal energy of a unit mass of 
the mixture. 



m a 



Since a is the number of molecules of combustible 
unit volume and Q,' the energy liberated in a chemical 



reaction of a single molecule of the combustible, aQ,' is 
the store of potential chemical energy of a unit volume 
of the mixture and aQ'/p is the amount of chemical energy 
of a unit mass of the mixture. Hence, equation (20) may 
be formulated thus: The sum of the thermal and chemical 
energies of a unit mass of the mixture during the com- 
bustion process remains constant. That this sum is the 
same before and after the combustion process seems obvious 
since it appears to be a simple consequence of the law of 
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If 



appears 
Thus , 



conservation of energy.. Since, however, we here have 
to deal with the phenomena of thermal: conductivity and 
diffusion equation (20), referring to different stages 
of the combustion process, in no wise 
of the law of conservation of energy 
actually the case), the. relation Dp 
coefficients of thermal conductivity and" diffusion were 
hot observed it would "be easy to show that e quation (20) 
would not "be correct. The sum of the energies at some 
points of the combustion zone would be greater than 



a consequence 
if, (as is 
between the 



c p *o + a 0 qi/p Q 



the 



and at others less. IPor this reason 
demonstration of equation (20) given by Zeldovich and 
Prank-Kamenet sky is an essential supplement of the work 
of Lewis who first postulated it. We shall see below how 
this lav/ can be generalized to the case where Bp is not 
equal to A/ Cp. 

Prom relation (20) Lewis did not derive any clear 
and simple conclusions for the general case of flame prop- 
agation. He computed, .with the aid of this equation, 
the flame propagation for the individual and very special 
case of combustion with the decomposition of ozone. Un- 
fortunately he evidently permitted important errors to 
enter in his derivation. It was possible for Zeldovich 
from this equation, together with all previous results, to 
arrive at general conclusions for the combustion theory 
which we here present, working them out to a greater 
accuracy. 



According to equation (20) 



a = 



f a-o _ c-p (T-To )1 



I" 



Po 



J 



and since 



t = 



IT 



- -^M 1 



- sir 



and 



hence 



a = 



HM 



qp_ Ti_ 



- T 



Ti - T f 



= a 



p T-t- T 



P o 1 - x 0 



= a, 



To 
T 



(21) 
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Let us return to formula (20) .giving the relation 
between the concentration of the c ombus t ijble a and the 
temperature in the heating and reaction zones, rewriting 
it in the following form: 



or 



am cp(Ti-T) _ , ^(yi) _ c p (Q? 1 -T)Mo _ 3?i ~ 

P o.t S, . " L " T i - 



m 



rn 0 
1 0 



v T, - T n s 



or 



B = ~~ = ( 1 - T) 



That is, the relative weight concentration of the combustible 
B at any point x of the combustion zone is equal to the 
difference between unity and the relative rise in temperature 
T at this point. Since. for x = B = 1 and T = 0 and 

for x = +<», B = 0 and T = 1 the curves 3 - x and T - x 
are symmetrical as schematically shown in figure 17. 

Having thus obtained the relation betvreen the number 
of fuel molecules and the temperature we may now find the 
velocity of flame propagation without making the incorrect 
assumption of the nondependence of the speed of reaction 
w on a which assumption was made in the derivation of 
formula (9). This computation was first carried out by 
Zeldovich and Frank-Kamenet sky . * 



We leave out of account the evidently incorrect com- 
putation of Lewis and von Elbe for the particular case of 
the flame propagation in ozone. 
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In the case of nonorn'olecular reactions 



-e/rt 



w = kae 



Substituting this expression in the integral I of formula 
(13) making use of (21) there is obtained 

2 m r'kae^ m _ E/ETi T r To e -ie/ET,« 
J a Q J T Ti - T 0 

The magnitude T 0 /T may be taken outside the integral sign 
and placed equal to T 0 /T 1 because T' is near in Talue 
to T 1 . In this case 



ke~ E / HT i T 



S\3 



T a - T 0 T x E 
Pi 

where j = / e~Pgdp = 1 - e~Pi(0i + 1) 
Cho os ing 



2HTx 3 3RT! 3 4ET 3 
0 = — — , , ■ 1 

1 S E E 

we obtain for p x the values 2, 3, 4, and for j the 
values 0.6, 0.8, 0.9. Thus, we dp not, make any large 
error in the expression for the velocity in setting j = 1. 
In this case, substituting the value of I in formula (14), 
we obtain , * 



These formulas differ from those of Zeldovich only in 
that under the square root' thVre st and' '"the magnitudes 
-T 0 /T 1 and (Tq/Tj) 8 for monomolecular and bimolecular 
reactions, respectively. The absence of these multipliers 
in the formula of Zeldovich is associated with a different 
manner of writing the equation of diffusion. 
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~ J ?o c V -Si (Ti - T 0 ) 2 V b / 



/feXcpk T 0 /RTj^ 3 \ 2 -E/RTj . , . 

_ / - i J I ±— ) e cm/sec (22) 

y p L3 Tj_ E 



The corresponding computation for the biaolecular reaction 
(w = ka s e -E/HT) g i veB 



M(kao)/To\ 2 e" E / RT i /HT^ V 



/ 4x(ka 0 )c- p a /ToY /B^V -E/ET- 



(23) 



Since the only magnitudes depending on the pressure 
at which the test is carried out are a 0 and p 0 ("both 
directly proportional to p) the linear velocity v 0 for 
monomolecular reactions is inversely proportional to ^p 
and for bimolecular reactions does not depend on the 
pressure. 

In the derivation of the above formulas we assumed 
that the main part of the reaction proceeds in a zone 
situated near the maximum temperature T x * and is limited 
to the temperature range 6 = T x - T 1 . As we have seen 

*The manner by which the combustion of all the initial gas 
supplied at a mass rate v m is maintained when the com- 
bustion proceeds in a mixture 90 percent diluted with the 
products of combustion and containing only 10-15 percent 
of initial number of molecules of the fuel may appear 
strange. As before this gas with low fuel content moves 
with mass velocity v m and it may be asked where the re- 
maining SO percent of the mixture is burned. There is no 
contradiction here and the explanation is that the remaining 
90 percent of the fuel is supplied to the boundary region 
not by the. motion of the gas but by 'the diffusion of the 
fuel molecules to the place of reaction. Thus all the fuel 
burns in the combustion zone but under conditions of strong 
dilution with the reaction products. 
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the method of solution is approximately true only for 
6 a ^ 0.25 T a .. Since the number of molecules a de- 
creases with increase in the temperature according to 
formula (21) the mixture, in the range of 6i (from 

- T 1 to Tj) will consist mainly of the reaction 
products and contain .a small percent of the initial 
regulating substances . In connection with this it is 
necessary to check to what extent, under these conditions, 
the assumption will be justified that the most rapid re- 
action takes place within the zone and to what extent 
the reaction may be neglected at temperatures less than 
T 1 = T j - 8 X while maintaining the conditi on 8 X 0.25 T x . 

For a bimolecular reaction the conditions will be 
least favorable. For this case 




where E is a certain constant depending on T and 
6 s Tj - I. If the reaction velocity is represented as 
a function of 6/T x the shape of the curve will depend 
on the value of y = B.T a /E. Carrying out the computation 
for y ~ 0.05, 0.1, and 0.2 we shall find that if the re- 
action velocity w may be neglected for y = 0.05, 
e/r x ^ 0.25, that is, consider that practically the entire 
reaction occurs at 8/Tj <> 0.25 then this is no longer 
true for 7 = 0.2. Hence for a bimolecular reaction the 
application of formula. (23) is restricted to the values 
RT3./E = 0.1 or values of the energy of activation E 
larger than 40,000 at maximum temperature T a S 2000°. 
Eor mono.molecular reactions conditions will be more 
favorable and formula (22) will be applicable to cases 
of practical interest. 

The above approximation, fundamental for the theory, 
may be analyzed more accurately though in a less clear 
fashion in the following manner: As we have seen the 
velocity of^flame propagation is determined by the in- 
tegral / w dT = I. We have already shown that this 
T • 
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integral should depend little on the choice of T 1 . In 
other words it should differ little from its •limiting 

value / w dT. In order that this he true it is nec- 
essary' that w(T') should be sufficiently small. We 
obtain a correct result only if the difference Q y = f 1 '- J S' 
thus determined does not exceed 0.25 T x . 




(24) 



Figure 18 shows graphically j as a function of for 
Y = 0.05,- 0.1 and 0.2. The arrows indicate the values of 
Pi for which j does not differ by more than 10 - 15 
percent of its limiting value. We see that these values of 
& 1 will be 4, 2.5, and 2 and the corresponding values 
of 6a/Ti s Ypi will be 0.2, 0.25, and 0.4. Thus, here 
too, we come to the same conclusion as before, namely, that 
the solution is valid only up to values of Y < 0.1 since 
for large values of Y the reaction zone is propagated in 
the temperature range B ^ > 0.25 T x and the assumption 
that the main part of the bimolecular reaction occurs in 
the temperature range near T x becomes untrue. We may 
note, by the way, that the limiting value of the integral 
«j = 2 obtained above (see formula (23)) is not correct and 
is' associated with neglecting the term Yg in comparison 
with unity in formula (24). The value of this integral 
depends in fact on the parameter Y and for the values, 
most often encountered in tests the limiting value of the 
integral is equal to unity. In other words under the square 
root sign in formula (23) there should stand 2 instead of 4. 

In the zone of reaction x > 0, as we have seen, 
Tj - T is small and of order of magnitude equal to RT^/E. 
We shall denote the concentration a of the fuel molecules 
for x > 0 by a eff . Then according to formula (21) a eff 
will be of order of magnitude equal to 
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Tp E a 0 c p /RT 1 a \ T 0 _ a 0 c p (.t RTj 8 /T 0 \ 

" a eff = ' a o « To l \ e / T i " -M 0 Q»N B \T i / 

Po c p RTg g fEp\ poCpTp RT i pc-p RTi s ' 



If we consider .that the fuel concentration in the 
combustion zone is constant and equal ,to a eff we are 
justified in the solution of the problem of velocity of 
flame propagation, in making use of formula (16) derived 
for reactions of zero order substituting S = ka e ff for 
monomolecular and ka 3 e ff for bimolecular reactions. We 
then, as may be shown without difficulty, obtain the same 
result as with the strict treatment expressed by formulas 
(22) and (23). Formula (16) and the expressions for the 
effective concentration &eff are thus sufficient for 
obtaining an expression for the velocity of flame propa- 
gation in various kinetic cases. In particular this re- 
mains true also for the case where the reaction occurs 
between two components a and b, for example, according 
to t he law . 

w = kabe 

If there is a deficiency in the component a and 
an excess in the component b the concentration b' of 
the component in the reaction zone may be assumed equal 
to its concentration in the. combustion products (it is 
computed by the stoichiometric equation from the initial 
state composition of the mixture) and the concentration 
a eff may be computed by formula (25). Then substituting 

the value w - ka ef f b « e~ E / RTl thus computed in formula 
(16) we obtain a correct result. It is interesting to 
note that the actual effective concentration of the fuel 
molecules in the reaction zone does not depend directly 
on the initial concentration a D but is determined only 
*y m ^^^f l ti t ^ en ?P e . r , a -* ure . of, combustion T x , . Thus, if 

the 'lit terms' maintained 'constant for various values of 
a 0 (by substituting, for example , one inert diluting 
substance for another with a different density or specific 
heat or by varying the initial temperature of the mixture) 
then the actual concentration of the fuel in the combustion 
zone remains constant. 
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Iff view of the above circumstance it would "be prac- 
tically impossible for us to determine the order of the 
reaction from flame propagation tests at constant pressure 
(varying a 0 by dilution with an inert gas). Since, how- 
ever, a e ff is proportional to the density of the gas or 
its pressure we may determine the order of the reaction 
from tests at various pressures. In the case of the com- 
bustion of a bimolecular system consisting of molecules 
A and B we may establish the order of the reaction only 
relative to that component which is in excess since only 
its concentration is entirely determined by the initial 
concentration and is equal to the difference between the 
initial concentration and the number, of molecules enter- 
ing into the combustion. Thus, in the example considered 
below of the oxidation of CO, from a comparison of the test 
data with the formulas for the propagation velocity we can 
assert only that the speed of reaction is proportional to 
the water vapor concentration, does not depend on 0 S in 
oxygen-rich mixtures and is proportional to the CO con- 
centration in mixtures with excess of the latter. More- 
over from, the dependence on the pressure it may be con- 
sidered that the reaction has a second order. If the 
assumption is made that the kinetic law is the same for 
mixtures rich in oxygen and rich in CO it follows from 
these data that the speed of reaction is proportional to 
the product (CO) (H s 0). If, however, this assumption is 
not made the following law would satisfy the test data: 
The reaction speed is proportional to the product of the 
general pressure by ( H s 0) in excess 0 3 and (CO) ( H 3 0) 
in excess CO. 

It should be stated that although the new theory, 
with the aid of the test results on the veloc it ies ■ of 
propagation, provides much information on the kinetics, 
of homogeneous reactions at temperatures fromwhich the 
reaction time is measured in 10~ 4 - 10 -5 seconds (pro- 
vided by no other methods ); from this peculiar nondepend- 
ence of the effective concentration on the initial con- 
centration; no complete guarantee of the uniqueness of the 
solution of the problem can always be given. 

In the derivation of f armulas (32) and (23). we made 
a number of restricting assumptions not always corresponding 
to the true conditions of the flame propagation. Our furthe* 
problem will be to free the theory of these restrictions 
and to extend the formulas obtained to any real case. 

These restricting assumptions were: 

(1) the specific heat was constant 
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.('2) . the rat io A / Cp , *Jhat is , the ratio of ' the 

heat conduct ivity ■ to the specif i'c ' heat | 
was constant 

(3) the number of molecules did not vary during 

the reaction 

(4) the coefficients of heat conduction and 

diffusion were equal 

In the case that V c p . is constant hut X is 
variable, it is necessary to employ as the variable not 
the temperature but the heat content H and therefore the 
equation of heat distribution in the form (!•)• The con- 
centration field a/p will be similar to the field of 
heat content, equation (21) being replaced by the equation 



Po H i - H c 



a = a o "f- P * ; (36) 



The method of solving the equation will be entirely anal- 
ogous to that which was earlier employed except that H 
will everywhere be substituted for T. Since, however, 

in the integral / w dH the temperature enters ex- 
^ i 

plicitly (w = Se~ E / R ®) and H ' is the heat content 
corresponding to a certain temperature T' near Ti it 
is necessary to transform this integral to the variable 

T« 

T in the form / wc p dT where Cp is variable. More- 

over, if the reaction is monomolecular w is proportional 
to a which is turn in linearly connected with H. There 
thus enters under the integral c p 3 and for bimolecular 
reactions c p 3 « However, since the temperature interval 
in which the reaction occurs is not large it may be assumed 
that within it Cp is constant corresponding to a tempera- 
ture near the combustion temperature and a composit ion -of 
the gas approximately that of the products' of combustion. 
We shall denote the value of Cp corresponding to the 
temperature T x by c* and the corresponding value of 
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the heat conductivity "by A*. Computation shows' that 
formulas (22) and (23) in their second form hold also 
for variable specific heat except that. for Cp and A 
th&re is substituted c* and A*, that is, the specific 
heat and heat conductivity of the products of reaction 
at the maximum combustion temperature T x . 

In the. case of non-constant A/cp and Dp we may, 
in integrating equation (11') assume A/c p = A*/c^, that 

is, consider it a constant magnitude because of the small- 
ness of the temperature interval T x - T'. VJith regard 
to the integration of equation (11) there is substituted 

'(cp/A) (T 1 - T 0 ) for (c p /A) ( T x - T 0 ) where (c p /Al is 

defined as 




Thus in the formula for the flame propagation in- 
stead of the magnitude Ac p for the monomolecular and 

Ac 3 for the biraolecular reaction there enter corre- 

P 

spondingly 



Due to the smallness of these corrections, however, we 
shall take no further account of them. 

If the reaction occurs with change in the number of 
molecules in the ratio ni/n 2 then in the reaction zone 
where the principal mass of the gas consists of the prod- 
ucts of reaction the number of fuel molecules will be not 

a a 0 Hi ■ - H 

~ = but 

P Po H l ~ H o 

a _ a 0 iij H 1 - H 

p " J~ Q ~ s Hl - E Q ' (2? 
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because the density of the gas varies in the ratio nj/n 3 . 
In correspondence with this '' n x /n 3 "and (n^/ng) 2 for mono- 
molecular and bimolecular reactions, respectively, enter 
the formula under the spare root sign. . 

Finally a very important case to consider is the 
one where ^/Cp * Dp their ratio being 



Dp = A : B 



This is particularly important in the cases of combustion 
of hydrogen in air or chlorine where the diffusion is 
determined "by the diffusion coefficient of hydrogen and 
the conductivity by that of the mixture. Simple computation 
shows that in this case it is necessary to introduce under 
the square root sign in the formula, the additional factor 
A/ B for monomol ecular and (A/B) 3 for bimolecular reactions. 

In the case X/ c^i Dp = A : B the fundamental equations 
(19) become 



= 0 




According to the law of conservation of energy in the 
initial and final components 

(<*>«, = a, = £2 = (0)00 = 8, = C P ^ ~ 
Po q« 

Since L = c p (T x - T Q ) = a 0 Q'/p 0 

In the reaction zone a and 9 , as before, approximately 
satisfy the equations 



B J£f + w = 0 A A!| +w = 0 
dx" 3 dx 8 
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Integrating in this interval we obtain 



00 




x" 



that is , 

^ = ^ and a S ; 6 + C 
dx B dx o 

Since according to the boundary conditions for 8 = 'B y) 
a = a 1 =6 1 we have 



B 

whence 



B 



In region I we may in equations (19) set w = 0 whence, 
integrating once, we obtain 



a £ - t„ e - o 

dx 



da 

B v m 

dx 



- v m a = 0 



Integrating a. second time we obtain 

a = a' e~ B and £ = 6 1 e A 
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where 6' and a' are the values of 8 and a for 
x = 0 "bounding regions I and II. 



a .o a ao 

— — - — = a; = a x 

P 0 P Po 



vrh^nce 



— = q, , - a = 8i - a 
P 



or by formula (38) 



- = ± ( 6, - 6 ) 

p B 

or 

A £l> , m m x A p (Ji - U *o - g f?q s 

a = P B ~ T) = B po a ° (Ti>T 0 ) = B S ° T V?; < 29 > 

In other words in. the region where the reaction proceeds, 
the values of a differ from those which would he obtained 
for A = B (see formula (21)) only by the factor A/B. 

In their final form the equations for the normal 
velocity of propagation, in the case of simple monomolecular 
or bimolecular reactions with one substance taking part 
(decomposition of C1 3 0, explosive substances), are written 
as follows: 

/ 2 A* c g k r T oY A Y n iV RT i s Y " sir , . 
v o = ^AiA-A— J e (30) 

for monomolecular reactions , and 
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for bimolecular reactions. 

It is necessary to say something about the widths of 
the heating and reaction zones. In the heating zone 
(x from -oo to 0) according to formula (7) for Tj - T 0 - T 2 - 

2 - T 0 = (T 1 - T Q ) e ax 



where a = v m Cp/A. This moans that the drop in temperature 
T - T 0 from its maximum value T 1 - T occurs e times 
in the distance 



b c P v m c pP 0 V o =V o Vo = V o T i = V o T o U ^ 



where v Q is the linear velocity of propagation and D 0 
the coefficient of diffusion at the initial temperature 
T Q . In the majority of cases at atmospheric pressure 

D 0 ~ 0.2; T 1 = 2000°; v Q = 50 cm/sec 

Hence the " heating zone 



£= "T^p x "^^ =2.6 X 10 2 cm=0.2-0.3 



mm 



The reaction zone | may be roughly estimated as follows 
from e qua. tion(ll) 



/dT\ = v^ 



where (dT/dx) 0 is the value of the derivative at the 
boundary of the reaction zone. We note further that for 
K = + 00 t. '(dT/dx)^ = 0. We shall consider that the- mean 
value of aT/dx in the reaction zone will be equal to 
s(dT/dx) 0 . Hence 
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— - i ZS££ /'• rjj " _ m V 
dx 2 \ v i ^o' 



dx . 2 X 
or, integrating we obtain 

As we see, the drop in temperature in the reaction zone 
is of the order (T x - T») = (T T - t g ) 0.2 whence the 
reaction zone is of the order of magnitude 



£' s 0.4 £ = 0,4 X 2.6 x 10" 3 £0.1 mm = 10" 2 



cm 



The mean free path is of the order 10" 5 cm at normal 
temperature and pressure. At 2000° it will be approxi- 
mately 10 times as large, that is, about 10~ 4 cm. This 
means that the entire width of the reaction zone is of 
the order of a hundred times the length of the mean free 
path. The number of impacts during diffusion in this zone 
is n = that is, of the order of tens of thousands. 

The time of stay in the reaction zone is 

T = ~— = — ■— |o = 10-JL =3>< 10 -» sec 
▼ v 0 Tj 50 20 

The time between impacts is l/u = 10~ 4 /l0 5 = 10~ 9 , that is, 
n = 10 4 . 

4. CONFIRMATION OP THE THEORY FOR SIMPLE REACTIONS 
(THEORY OP COMBUSTION OP EXPLOSIVE SUBSTANCES 
AND IN PARTICULAR N ITROG-LYC OL) 

Belayev (reference 4) proposed a theory according to 
which all easily vaporizable explosive substances' burn in 
the gaseous phase. The surface' of the explosive burns con- 
tinuously by the heat of the combustion conducted to the 
surface, the heat bringing about vaporization and raising 
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the temperature at the surface to the "boiling temperature 
J-j, and thus supplying a continuous flow of vapor to the 
combustion zone. For all combustible explosive substances 
capable of being vaporized vrithout igniting .( and which 
not only detonate but burn) this new point of view regarding 
their combustion appears to be much more probable than 
the old conceptions of direct combustion in the condensed 
phase. Belayev proved his hypothesis by direct test. By 
photographing the flame of nitroglycol he showed that 
between meniscus of the liquid and the flame zone there is 
a narrovr dark zone, the width of which is measured in tenths 
of millimeters at atmospheric pressure and attains . 1 milli- 
meter at atmospheric pressure and attains 1 millimeter for 
combustion at loitered pressure. It is thus possible to con- 
sider the theory of Belayev as proven by direct experiment. 
The mass rate of combustion v m may be directly measured 
from the rate of lowering of the meniscus (u) in the com- 
bustion of nitroglycol. This directly measured magnitude 
may be provisionally denoted as the linear velocity of com- 
bustion of the liquid nitroglycol. It would be more correct 
to call it the linear velocity of vaporization in the com- 
bustion of nitroglycol. The mass rate of vaporization v m 
is evidently equal to the mass rate of combustion of the 
nitroglycol vapor since all the nitroglycol burns in the 
form, of vapor. The linear speed of combustion of the 
nitroglycol vapor is therefore 



0 Po p n 



where (p) is the density of the liquid nitroglycol, p 0 
the density of the nitroglycol vapor at the test pressure 
and the initial temperature. 

Belayev applied the theory of Zeldovi ch-3Tr ank-Kamene t- 
sky to the computation ,of the speed of combustion of the 
nitroglycol vapor. We shall present the computation, making 
use of the equation of heat propagation, in the form of 
equation (2) . 

As was done before we shall divide, the entire process 
into that of region II where the main part of the reaction 
takes place and where it is possible to neglect the term 
v m Cp dT/dx and region I or the preheating region where it 
is possible to neglect the speed of reaction w. The 
difference between the present case and that described 
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above is that region I itself roust now be further sub- 
divided into region I « (the dark preheating zone- of the 
vapor) and region I»« of the preheating of the liquid 
explosive from the initial temperature T Q to the boil- 
ing temperature at the boundary of the liquid and vapor 
phase (fig-. 19). In both regions the same equation 
applies 

\ d 3 T ' dT 

x apr - v p — = o 

except that the constants X and c p are different. For 
the liquid these magnitudes will be denoted (X) and (c p ). 
Since the vaporization associated with the expenditure of 
thermal energy occurs at the boundary between the- liquid 
and vapor phases of the explosive subs t anc e, the flow of 
heat. X(dt/dx) t » from the vapor side of the boundary sur- 
face will be greater than the flow (A) (dl/dx)^' I from 

the other side of the boundEii-y surface in the liquid phase 
by the amount of energy required in vaporizing v m grams 
of the substance per second per unit area of surface. 
Hence, 



where P is the latent heat of vaporization of one gram 
o f t h e e xp 1 o s i v e . 

The temperature distribution is qualitatively in- 
dicated in figure 19. Integrating the equation for region 
I ' we obtain 



Integrating the equation for region I*' \ie obtain 
M SDV = V -^p] (Tboil " *o> 
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whence for region I* we obtain 



^ (dl)' = Vm ( Cp iT ~ T ^ oil) + C C p] (T D0il - T 0 ) + ?] (34) 



or integrating, the width of the dark preheating zone is 
obtained as 

t - X ~i C P (g i I T boil) + [ c p3^boil ~ *o)+* 



The magnitude in the numerator under the logarithm sign 
is no other than the heat of combustion of 1 gram of the 
liquid explosive. Thus 



As regards flame propagation velocity the latter 
will be expressed by the same formulas (50) and (31) 
as for the gas because in zone II the equation is the 
same as for the gaseous case and the expression for the 
derivative according to formula (34) is 



dT y 

•dx J Q 



A study of the speed of combustion of nitroglycol 
shows that (u) increases in direct proportion to the 
external pressure p, that is, (u) = p v /(p) where ' 
p 0 is directly proportional to p ana the density 
of the liquid (p) does not depend on the pressure. As 
we have seen, the . formula v 0 varies in inverse proportion 
to J p for monomolecular reactions and is independent of 
p for bimolecular reactions. Thus, in order to satisfy 
the test data with regard to the dependence of the speed 
of combustion of liquid nitroglycol on the pressure, we 
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must assume that the reaction is bimolecular . On the 
other hand, direct study by Appin of the slow decomposition 
of nitroglycol vapor at low temperatures and pressures 
has shown that the spood of reaction obeys the monomolec- 

35000 



ular law w = 10 14 ae - 



KT 



Belayev assumes that at the high temperatures correc 
spending to the combustion zone, the speed of activation 
of the molecules as a result of the impacts becomes in- 
sufficient for the establishing of a stationary concen- 
tration (computed from statistical data) of the active 
molecules. The speed of reacti-on then "begins to he deter- 
mined by the number of activating impacts, that is, becomes 
bimolecular. A similar phenomenon has long been known for 
all monomolecular reactions which with the lowering of 
pressure show a drop of the molecular constant and hence 
for sufficiently low pressures become bimolecular. It is 
also known that with increase in the temperature the drop 
in the constant with decreasing pressure becomes more sharp 
and thus the pressures at which the reaction passes over 
into the bimilecular increase. 

It is possible that in this way the hypothesis of 
Belayev should be explained according to which at very 
high temperatures of the reaction the decomposition of 
the nitroglycol is bimolecular, even at atmospheric 
pressure. It should be noted that according to test 
results on the kinetics of monomolecular decomposition 
the reaction products as a rule are capable of as'good : 
activation of the molecules of the initial substance on 
impact as the impingement of two molecules of the initial 
substance. Thus the bimolecular ity of these reactions is 
exceptional. Whereas, for example, for the usual bimolec- 
ular reaction of the type 2HJ = H 3 + J 3f the decomposition 
is possible only on the impact of two particles of the 
initial substance and w = ka s e" 3 / SI ; in this case the 
activated molecule of the nitroglycol decomposes by itself 
and the reaction is possible on impact iirith any particle, 
particularly with particles of the products of its de- 
composition, whence 



w = Z a(M)e- E / Rn? 

where (M) is the number of all molecules (both of nitro- 
glycol and of the products of its decomposition)' in a unit 
volume. 
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The number (M) is determined by the initial number 
of molecule's of nitroglycol a 0 and is equal to a Q T 0 /!r 1 ' 
whence 

w = (Za 0 ) ae" E / RT ^ 

Remembering the character of the derivation of the 
expression for the velocity v Q we shall see that it 
will be determined not by the expression for bimolecular 
reactions but by the expression for monomolecular reac- 
tions with this difference only, that instead of the 
constant of first order k( ~ 10 13 - 10 14 ),we shall have 



Za = y2T7(T 2 ua n 1?- = 10 » 
In the case that A = B, we obtain 



° y po La w e / 1 (36) 

As in the case of purely bimolecular reactions v 0 
doss not depend on the pressure p because a 0 and p Q 
"are directly proportional to the pressure. Hence u and 
v m will in correspondence with the test be directly 
proportional to the pressure. Since p 0 = p,a 0 /M where 
p, is the molecular weight of the nitroglycol, 



The fundamental magnitude for the computation of the 
speed of combustion is the heat of decomposition L of the 
nitroglycol in the flame during its combustion in an inert 
atmosphere. In the case of its combustion in air there 
arises, as we have seen, a secondary flame of combustion of 
the products of the primary decomposition. This secondary 
process has no relation %tfith the speed of combustion and 
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we are therefore interested only in the heat of combustion 
of the primary process which readily occurs- in an • inert 
atmosphere. Direct measurement of L has not teen, made, 
Appin, however, determined the products of the primary 
flame "by direct analysis. The latter gave the following 
final chemical equation for the decomposition: 

0 3 H 4 (ON0 3 ) a = 2N0 + 1.7C0 + 1.7H 3 0 + Q.3C0 3 + 0.3H 3 

that the greater part of the molecules is 
according to the equation 

C 3 H 4 (03r0 3 ) B = 2110 + 2C0 + 2H 3 0 

and a smaller part according to the equation 

C 3 H 4 (0N0 3 ) s = 2H0 ■+ . 200 + 2H S 

Computing the heat of this decomposition of nitroglycol 
on the basis of the heat of formation, as recommended 
by A, Schrait , Belayev arrives at a value L = 450 calories 
per gram. Knowing L and using the mean specific heat 
of the reaction prodxicts given in the tables of Lev/is 
and von Elbe, Belayev arrives at the following value for 
the maximum combustion temperature: 

Tj - T 0 = 1350° or 0?! = 1350°+ 300° = 1650°K 

The energy of activation E is taken by Belayev from the 
tests of Appin for low^-temparature decomposition of the 
nitroglycol vapor equal to 35,000 calories. The specific 
heat of the reaction products for Ti = 1650°K is computed 
by Belayev from the data of Lev/is and von Elbe as c|/ =0.35. 
The heat conductivity of the products for T = T x was com- 
puted from the mean conductivities at 0°C and then corrected 
for the temperature 1650°K by the correction of Ceteriand, 
The ratio of the number of initial to the number of final 
molecules n x /n 3 = 1/6 (according to the stoichiometric 
equation). 



This shows 
dec ompos ed 
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We shall carry out these computations "by a 'somewhat 
different method. According to the kinetic theory of 

gases — si ^ 3c ff (for derivation see below, section 5).. 
Po 3 n o 

In the case of the decomposition of the vapor of pure 
nitroglycol n 0 = a 0 

u 3 ~ 3 X 83 X 10 6 X T x ~ 3 X 83 X 10 6 X 1650 fi 4 X 10 11 
p. jj. M- 

Both as regards the number' of impacts and the heat 
conduction the principal role in the combustion zone will 
be played by the molecules of the decomposition products 
of nitroglycol, these being the lighter. Their value of 
ji must therefore enter the formula. Assuming = 30 
we have 

u 2 = | x 10 10 ; c p = 0.35 



whence 

M _ 1.5 X 10 9 
Po a o 



E ' 35000 

; e" RT, = . e 3300 = ,-io.e = 2 . 5x 1Q -5 



fl° Y - (_300\ a „i.£i. I (Bl\ S - f 2:x 165Q3V = 2 34 x 10 4 
KtJ ~ Vl65 0 ; SO'^-bU "V 35000 J" 2 ' 34 * 10 

Substituting all these values in formula (36) we obtain 



k x 1. 5 x 10 9 2. 34 x 10 4 X 2. 5 x 10~ 5 ~ _ , 

v = / s — ■ = /50 = 7 cm/sec 

0/2 x 10 5 30 x 6 v ' 



and the corresponding mass rate of reaction v m = v Q p Q 

at p = 7 x 10~ 3 gives v m = 4 x 9 x 10~ 3 g/sec cms. The value 
actually observed, by Belayev at room temperature is 
4 x 5 x 10~" 3 g/sec cm2. The agreement between, the theory 
and experiment is therefore extremely good. 
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On increasing the initial temperature T Q at a 
given heat" L the maximum temperature Tj increase's and 
hence also the velocity of reaction and the velocitjr of 
flame propagation. Belayev computed temperatures Ti 
corresponding to various values of T Q from 20° to 200°. 
Moreover, he experimentally determined the speed of com- 
bustion u (proportional to v m ) at various temperatures 
T 0 . We present below some of the results. 



TABLE 2 



T o .in °C <S 1 in °K T in °C u = Zm 

Po 

20 1650 20 0.29 

50 1685 60 .33 

110 1753 100 .44 

140 1790 140 .52 

170 1828 180 .57 

200 1867 



The relation "between T 1 and u = v m/po ' as is 
seen from, the formula, is " ' ' 



ln ^ = ilk + C * ^ lg t: = 2X2T X + C * 



Substituting for each given T Q the computed T x and 
the experimentally determined u, we should obtain a 
linear relation between lg u/ T x and l/T a with the 
s 1 op e . 

^3 = ^1 - «00 



Figure 20 shows the experimentally determined relation 
between lg u/Tj and 1/^ for liquid nitroglycol (data 
of Belayev). We see that the points with some scatter 



*T 0 does not enter since, according to formula (36) with 
increase in T Q , there is a simultaneous decrease in a 0 
so that (a 0 T 0 ) s remains constant. 
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arrange themselves on a straight line with a slope equal 
to 4100. 

The above results thus also confirm the theory of 
flame propagation and give, for the energy of activation, 
a value practically coinciding with the value 35,000 
assumed bjr us according to the tests of Appin for the 
computation of the absolute velocities. 

In conclusion we shall compare the observed and com- 
puted widths of the heating zone. Substituting in formula 
(35) the value A* » 2 x 10 -4 . c| - 0.4 and v m = 4.5 x 10-2 
the heat of vaporization of nitroglycol F = 43,000 calories 
and the boiling temperature T-5 = 200°C. we obtain 
'£ =0.2 millimeters. The width of the dark zone at atmos-. 
pheric pressure directly measured by Belayev (microphoto- 
graph and flame photograph) gave a value of the same order 
of magnitude. 

Since the velocity v m varies in direct proportion 
with the pressure, with decrease in the latter, the width 
of the dark zone should increase in inverse proportion 
to the pressure and, as a matter of fact, the tests of 
Belayev for pressures of 100 millimeters give a w: dth of 
the dark zone of the order of 1 millimeter. Thus the tests 
of Belayev give excellent confirmation of the new theory 
of flame propagation. 

5. CHECK OF THE THEORY BY EXAMPLE OF A CHAIN OXIDIZING 
REACTION (THEORY OF FLAME PROPAGATION IN CARBON MONOXIDE , 
AIR, AND CARBON MONOXIDE - OXYGEN MIXTURES) 

One of the most important combustion reactions is 
that of the combustion of carbon monoxide. Together with 
the oxidation of hydrogen this reaction determines also the 
combustion of hydrocarbons (since it may be assumed that 
in the heating zone the hydrocarbons are first transformed 
into CO and H 3 which then burn) . 

The oxidizing reaction of hydrogen represents a chain re- 
action with many branches. (The question as to the applica- 
tion of the theory of Zeldovich to this. type of reaction still 
remains open.) In contrast to this the oxidizing reaction 
of CO represents a chain reaction, the velocity of which 
very soon attains stationary values, and we may assume that 
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notwithstanding the small width, of the flame front the 
velocity ofTfiis'f eaotion is completely determined by 
the temperature and the concentration at any given place 
within the combustion zone. This behavior of the oxi- 
dizing reaction of CO is associated with the fact that 
it takes place only in the presence of vtrater vapor, the 
products of whose decomposition (OH, H) are active cen- 
ters. Since jthe quantity of OH and H its limited by the 
moisture content the reaction velocity cannot exceed a 
certain definite value. In section 1 we enumerated, in 
dbts.il, the various properties of the flame propagation 
velocity in CO mixtures. Zeldovich analyzed the data 
of various foreign authors and also those of our coworker, 
Barsky, and showed that the properties of the CO flame 
and its absolute velocity may be obtained from our general 
theoretical formulas if it is assumed that the oxidizing 
reaction velocity of -CO is determined by the following 
kinetic law: 



a[oo a 3 i[° 3 ] a [bo] 

TZ = - 2~TT = - — T7 — 



dt 



dt 



dt 



Ze 



3 5000 
ET 



[H 3 0][C0] 



(37) 



We may note that H 3 0 appears only as a catalyzer 
and is not required in the course of the reaction. In 
order to present more clearly the reaction velocity 
according to (37) we construct a table of values of 

-25000/RT , 

(JU = e (the probability of the reaction computed 

for one impact of the molecules H s 0 and CO) and the 
reaction times T (time during which the products 0 S 
and CO are decreased to l/e of their initial values).. 



TABLE 3 



T° K ' 


U> X 10 3 


T 


sec 


1300 


0. 03 


3 


x 10~ 3 


1600 


. 33 


3 


X 10" 4 ' 


2000 


2. 8 


3. 5 


x 10~ 5 


2400 


7 


1.4 


x 10~ 5 


3000' 


'" 25 




x 10- 6 



In order to compute the velocity of flame propa- 
gation according t o ' the kinetic law of equation (37) 
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we shall distinguish, case 1 when the initial gas contains 
00 in considerable excess of the stoichiometric quantity 
and case 2 when, there is a deficiency of CO. 

Case 1.- As we hare already seen above, the reaction 
in. the combustion zone occurs in a gas, the composition of 
which is. near that of the final products. In the cooled 
final products the concentration of carbon monoxide will 
be (C0) a = (C0) Q - (C0 2 ) 1 where (C0) Q is the number 
of molecules of CO in a unit .volume of the initial mix- 
ture and (COg^ is the number of molecules of C0 3 in 
a unit of volume of the final products (after their cool- 
ing to temperature T 0 ) or, what amounts to the same thing, 
the number of burned molecules of CO. In the combustion 
zone where the temperature is near the maximum temperature 
of combustion 3? • the number of CO molecules in unit' 
volume (the actual concentration of CO) will be 



[co] = [co], 

The number of molecules of H 3 0 in the combustion zone 
varies in relation to (H 3 0) 0 only in the ratio Tq/Tj 
since H 3 0 is not required in the reaction. Hence the 
actual concentration of H 3 0 will be 



[E 3 0] = ±3. [H3 0] 0 



Thus 

„ = Ze - E / R0? (^) 2 CH 3 0] o [C0] (38) 

We soo that in the above expression the oxygen concen- 
tration, that is, the only concentration which varies 
considerably in the combustion zone does not enter. Thus 
in realtion t.o 0 3 the reaction is of the first order and 
therefore the mathematical expression for the velocity 
of the flame propagation will be similar to formula (16). 



For the sake of clearness, however, we shall derive 
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/ 



v o ~ ~T J 2X( ^' / w d * = ~- /2Aq' / w dT (39) 



Po* V J ~ Po L 

nvi 0 



By substituting w from formula (38) and carrying out 
the corresponding computations, we find 



/ 

v 0 = / 1400 — — - e 1 * 19 V a ' " (40) 

where I = 25, 000. 

We shall derive formula (40). -Substituting w from 
formula (38) under the integral sign of formula (39) we 
obta,in 

W dT " 2 "If" \tj ° CH 2 0] o [CO], (41) 



- to xz <3lL 1 ETi 3 /T 0 \ a 



According to the- kinetic theory of gases the heat con* 
ductivity is 

A = 3 l u P u 



where c is the molecular specific heat of the mixture 
(at the combustion temperature Q ~9) 
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Z = /ina'u .-ft 



whore n is the number of molecules in a unit volume of 
the mixture at the combustion temperature and equal to 
n 0 T 0 /3?i (where n 0 is the number of molecules in the 
cold mixture). Por T = 0^ the thermal velocity of the 
molecules is 



u = u 



300 



where u 300 
t empe ratur e , 



is the velocity of the molecules at room 



that is, for our mixture 



u 



300 



= 5 x 10 cm/sec. 



AZ 1 u g C 0?! p ■ 1 u 3 C 
Po 3 V> n o T o Po 3 M- n o 



75 x 10 s 

300 nn 0 



L Po 



[CP a ] 1 



since ^'(CO;^ = Lp 0 is the heat 
of a unit volume of the mixture. 



given up in the combustion 



L = 



where M Q is the ratio of the weight of the burning CO 
molecules to the weight of all molecules in the initial 
mixture, that is, in view of the nearness of the molecular 



weights of 
is equal to 



CO, 0 2 
(COs),. 



and N 



3 



and because .the burned 



w 0 - 5 

0 



and hence 



l qTcoTJT 
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where the heat of formation of a molecule of 0 0 3 is Q = 
7 x 10 4 -calorics. Substituting all the. obtained values 
in formula (42) we obtain formula (40) taking a mean value 
of {j. Z 30. 

Ras e 2.~ For a deficiency of CO its concentration 
in the ccmbustion zone is small and rapidly changes along 
the zone. Since the reaction velocity is proportional to 
(CO) we should here make use of the same derivation that 
was made above for monomolecular reactions. According to 
formula (21) the CO concentration in the combustion zone 
is 



[CO] = [C0] o 



Substituting this expression in formula (37) and inte- 
grating we readily obtain 



\ dT Z Sg- (|*)' Ze" E/llTl [00] o [H 3 0] o (43) 



•7 1 



Comparing this expression with (41) we see that the former 
differs only in that for the concentration (CO), there is 

C HI 1 ! * 

substituted the magnitude (C0) eff = ' (C0) Q and 

(C0) 0 p 

since for deficiency of CO ■ = — ^5 

L q' 



ou ef f i e 

3for monomolecurlar reactions, however , as we have seen, 
this magnitude being multiplied by T 0 /T a corresponds to 
the actual concentration (CO) in the reaction zone. Thus , 
the f»'=!x!la for the flame propagation velocity in mixtures 
v-ich a deficiency of CO is the same as for the case of a 
considerable excess of CO except that instead of (CO) 
there is substituted (C0) eff . Practically, however, a 
considerable difference is obtained, Whereas in rich 
CO mixtures' (CO^ increases with increase in (C0) o in 
rich 0 3 mixtures (C0) eff as is seen from formula (44) 
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does not depend on (C0) o . if the combustion temperature 
is maintained constant (for example, with the aid of heat- 
ing of the mixture by some method).. 

We thus arrive at the conclusion that the velocity 
of flame propagation in mixtures with equal combustion 
temperature remains constant for all mixtures with excess 
0 3 and is independent of the excess of 0 3 . This con- 
clusion was confirmed by our coworker, Par sky, who measured 
the velocity of flame propagation in mixtures with various 
excesses of oxygen and the same concentration (CO) (vary- 
ing the ratio 0 3 /N 3 in the mixture and thus maintaining 
a constant combustion temperature). The flame propagation 
velocity actually did not depend on the excess oxygen. On 
varying by the same method the CO content in mixtures 
with excess CO and keeping the quantity of oxygen the 
same (that is, again maintaining Ti constant), Barsky, 
in correspondence with formula (40), showed that the veloc- 
ity of flame propagation is pr opo rti onal to the square root 
of the excess CO, that is, J (C0) 1 

Jot a simult riation of both components the 

combustion temperature must necessarily likewise vary and 
it is •impossible by such simple method to bring out the 
dependence of the reaction velocity on CO and 0 3 . Ex- 
periments of this kind, conducted by other authors, are 
fundamentally in error. By utilizing, however, a certain 
in i t ial . preheating of the mixture it is possible to maintain 
T x constant on varying the ratio of the components. Making 
use of the data of Passauer and others (reference 5) Zel- 
dovich conducted such a test of the theory and arrived at 
the same result, namely, that the propagation velocity at 
constant T, does not d epend on 0 3 for an excess of 0 3 
and is proportional to J ( C 0) a for excess CO. According 
to formula (40) the propagation velocity is proportional 

to J (H 3 0) q (that is, to the square root of the concen- 
tration or the partial pressure of the water vapor). 

On replotting the results of Fi ock an d Marvin, shown 
in figure 15, in coordinates v - J (H 3 0) 0 we obtain good 
straight lines in agreement with the theory (figure 21). 
Por large water vapor contents there are deviations for 
the reason at least that the thermal character of the 
mixture then begins to change., 

As regards the dependence of the propagation velocity 
on the pressure of the mixture it follows from formula (40) 
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that for constant percent of moisture content the prop- 
agation velocity does not depend on the pressure -(because 
(CO^, (C0 S ), (H 3 0) 0 are proportional to p) and at con- 
stant partial vapor pressure v Q is inversely proportional 
to v ^p. ¥ e know that just this kind of dependence on the 
pressure is observed in experiment. 

It was most interesting to compare the' computed 
absolute value of v 0 for mixtures of various ratios 
of CO to air. The computation led to the results shown 
by the dotted curve of figure 10. There are, to be sure, 
certain deviations from the test results' in the region of 
the maximum but these, however, are of the same order of 
magnitude as the deviations in the test results. A similar 
computation for mixtures of CO wi.th nitrogen and oxygen 
in various proportions likewise led to good agreement with 
the test results of Jahn given in figure 12. In general, 
it may be assumed said, that the approximate theory gives 
a fair description of the test results in this case, 

A comparison was furthermore made between the com- 
puted and observed propagation velocities for various 
initial temperatures T Q for various CO air mixtures. 
It was found that the deviations between theory and ex- 
periment in this case are not large. A qualitative ex- 
planation was also obtained of the fact that the increase 
of the propagation velocity with the pressure in oxygen 
mixtures is small. The fact is that whereas it is per- 
missible with air mixtures not to take into account the 
incomplete comtustion -associated with dissociation, in 
the case of oxygen mixtures where temperatures up to 
3000°X and above are developed, it is necessary, to take 
dissociation into account. Details of the computation 
will be found in our article (reference 6) while here 
we give the result only. The computation leads to an 
increase of the velocity with the pressure as p n where 
n = 0.12; whereas experiment gives the increase of the 
velocity with pressure as p Q where m s 0.2. Thus in 
this case (if the tests are considered reliable) the 
theory gives only the qualitatively correct result, namely, 
a small increase of v Q with the pressure in rich oxygen 
mixtures. 

Summarizing all that has been said above we may state 
that the new theory of combustion, through the choice' of the 
form of kinetic law for the rate of oxidation of CO, which 
law contains only a single arbitrary constant (E = 25,000), 
rationally accounts for all the numerous experimental re- 
sults on the flame propagation of CO in air and in oxygen, 
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quantitatively in most cases and qualitatively in only 
a few cases. It is very int erest ing . to confirm directly 
the correctness of the kinetic law postulated by us for 
the r ate of oxidation, of CO. Ho experimental method, 
unfortunately, has up to the present time "been developed 
for the study of the reaction kinetics at high temperatures 
such that the reaction proceeds halfway in a time interval 
of 10~ 3 - 10~ 4 seconds. We shall, therefore, proceed by 
a logical method and by examining closely the question of 
the mechanism of the oxidation of CO, attempt to find 
the required kinetic law. 

It has been assumed that the catalytic action of 
water vapor in the oxidation of CO is based on a con- 
version reaction of the water vapor: H 3 0 + CO = C0 3 + H 3 . 
The hydrogen thus obtained is soon oxidized and therefore 
the total oxidizing reaction proceeds with the reaction 
velocity of a conversion reaction. Assuming that the 
conversion proceeds according to the simple bimolecular 
law we very simply arrive at the required law for the 
oxidation of CO: 

w = Ze,- E / RT CC0][H 3 0] 

Assuming the activation energy of the conversion reaction 
equal to 25, COO calories we obtain the law (37) employed 
by us. Unfortunately the velocity of conversion computed 
by the above formula does not at all correspond to the 
meager experimental facts which are available for the 
study of a homogeneous conversion reaction. Such a re- 
action in technical procedure,- as is known, occurs \*ith 
catalyzers and therefore few have concerned themselves 
with homogeneous conversion. On the velocity of conversion 
there are available data by Thompson (reference 7) and 
Kondratyev (reference 8) for reactions in quartz vessels 
at T = 1000°K. For the conversion velocity under these 
conditions Thompson obtains e ~ E /-^ = 2.5 x 10 -14 and 
Kondratyev a value one hundred times as large, namely, 
10~ 12 . According to formula. (37) the value obtained 
is approximately 10~ 6 , that is, 10 8 times as large 
compared with the results of Thompson and 10 6 times as 
large compared with the results of Kondratyev. The re- 
sults of Thompson should be considered as more correct 
because he experimented with a vessel of greater radius 
and the possibility of catalysis by the wall was less 
probable than was the case with Kondratyev. In any case 



'cans* 



KACA Technical Memorandum STo* 1026 71 



wo sec there can he no question of such a simple ex- 
planation of the lav;- (37). ..... ... 

Investigations of the velocity of the oxidation 
reaction, of CO at temperatures of the order of 1000 K 
in the presence of. water vapor show that this reaction 
is complicated by the chain process, the velocity being 
proportional to (C0)(H 2 0) and inversely proportional to 
{0g)» It is very probable, however, that under these 
conditions the primary active centers required for the 
start of the chain arise at the walls of the vessel as 
a result of a heterogeneous process, which, of course, does 
not occtir at high temperatures and. especially in flames. 
For this reason the extrapolation of the data obtained 
under the above-mentioned conditions to the flame will 
hardly he justified. Moreover, reactions in closed 
vessels have received but little study. We are fully 
justified in stating, however, that a reaction which 
appears as a chain reaction at T s 1000° remains a chain 
reaction at any high temperature whatever. It is well 
known, moreover, that the flames of moist CO reveal 
in the spectrum the OH radical. Kondratyev further- 
more, quantitatively measured (by the method of light 
absorbtion) the concentration of the OH radicals in 
rarified flames in moist CO having a temperature of 
1000° - 1200° K and showed that the concentration of OK 
exceeds its equilibrium value a thousand times. These 
facts clearly indicate that the OH radical is one of 
the -active centers of the chain oxidizing reaction of CO. 

Analyzing all possible chain mechanisms of this type, 
wc arrived at the conclusion that only one of them is 
logically possible at high temperatures. This chain con- 
sists of the following steps of elementary reactions: 



1. H + 0 S a OH + 0 N 

( Chain 

2. CO + OH = C0 3 + H ■> 

3. 0 + CO = C0 2 

4. 0 + H s - = OH + H Branched chains 

5. OH + H 3 = H 2 0 + H 



As a result of the branching of the chain the total 
oxidizing reaction will be accelerated in time until the 
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concentrations of OH and H reach vala.es such that the 
reverse reaction of (4) OH + H-*"-?>H 3 + 0 balances the 
direct reaction 0 + H 3 -5>0H + H. After this a stationary- 
oxidizing reaction is established. It is easy to show that 
the direct reactions 4 and 5 will be in equilibrium with 
the reverse reactions, that is, in the reaction zone there 
will occur the equilibriums 

4. 0 ■+ H 2 <=^0H + H 

and 

. ' 5. OH + H s <-==^H 3 0 + H 

The equilibrium constants K 4 and K 5 of these reactions 
are well known. + 

As regards the constants Kj , K 3 , and Z 3 corre- 
sponding to the velocities of reaction 1, 2, and 3 the 
situation is less favorable. For the constants Kx and 
K 3 there are approximately known the activation energies 
E a = 26,000 * 4000, E 3 = 10,000 ± 3000. The coefficient 
before the exponential for the first reaction is not known 
and for the second is of the order 0.01Z. ITo study was 
made at all of the constant K 3 . If it is assumed that 
reaction 3 occurs for every third impact then at atmospheric 
pressure Z 3 £ 10~ 3 . If the most probable values of the 
constants from the- experimental and theoretical points of 
view are taken, the series of reactions leads to a kinetic 
lav/ which agrees with expression (37) in the sense of the 
dependence of the velocity on the temperature and on the 
concentrations (CO) and (H 3 0) ; (it gives a practical non- 
dependence on (0 3 ))» For mixtures with excess CO, however, 
it gives an absolute reaction velocity about 10 times less 
than the actual. Since the measured values of the kinetic 
constants of the elementary reactions are bad the matter 
is evidently associated with the not entirely correct choice 
of these constants. In general, we may say, that the check 
of the new theory on the oxidation of CO leads to satis- 
factory results. 



Translation by S. 
National Advisory 
for Aeronautics. 



Hei s s , 
Commi tt ee 
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Continuation of figures 11, 12, 15 



Key for Figure 11 



1 - 87# 3ff 3 +l3# 0 S ; 

3 - 79$ 1 3 +21$ 0 S ; 

5 - 60$ U 3 +4C$ 0 3 ; 

7 - 20.J6 N 3 +80$ 0 3 ; 



2 - 83$ N 3 + 17$ 0 3 ; 

4 - 70$ 1T 3 + 30$ 0 3 ; 

6 - 4056 JT 3 + 60$ 0 3 ; 

8 - 1.5J6 N 3 +98.5$ 0 3 ; 



Key for Figure 12 



1 - 87.5$ K 3 +12.5$ 0 3 ; 

3 - 82.5$ 3!J 3 +17.5$ 0 3 ; 

5 - 75$ N 3 +25$ 0 3 ; 

7 - 65$ ar 3 +S5$ o s ; 

9 - 50$ U 3 +50$ 0 3 ; 
11 - 30$ K 3 +70$ 0 3 ; 
13 - 10$ U 3 +90$ 0^ ; 



2 - 85$ N 3 +15$ 0 3 
4 - 79$ N 3 +2i$ 0 3 
6 - 70$ N 3 +30$ 0 3 
8 - 60$ 3ff 3 +40$ 0 3 
10 - 40$ U 2 +60$ 0 3 
12 - 20$ N3+8O56 0 3 
14 -1.5 N 3 +98.5$ 0 3 



Key for Figure 15 



1 


- PC0+0 3 




100 


mm ; 


2 


~ Pco+o 3 




150 


mm ; 


3 


- Pco+o 3 




200 


mm ; 


4 


~ Pco+o 3 




300 


mm ; 


5 


** £co+o 3 




7 60 


mm ; 



NACA Technical Memoranda Ho. 1026 






Tigs. 1,7,8,9,10,11 
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Figure 1. 



Figure 7. 



Figure 8. Figure 9. 
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Figure 10.- Normal velocity of 
flame propagation 
of carton monoxide-air mixtures. 

Full curve - tests of Passauer, 
x - tests of Chitrin, 
o - tests of Baraky, 
o - tests of Jahn. The dotted 
curve was computed "by the 
formulas of the new theory of 
the flame propagation with 
the kinetic law (37). 



Figure 11.- Normal velocity of 

flame propagation in 
mixtures of CO with atmospheric 
Ng +- Og. The percent composition 
of No + 0 2 of each curve is 
constant (Jahn) . 
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Figs. 3,3,4,5,6 




Figure 3.- Wave front moving 

in a horizontal 
tube (d = 5 cm). The burning 
mixture consists of 55 
percent 00 and 45 percent 
air. At point a flame 
oscillations arise, 
(photograph by Barsky) . 




Figure 5 





; 







Figure 6 
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Jigs. 12,13,14,15 




Figure 12.- Normal velocity of flame 

propagation in mixtures 
of hydrogen with atmospheric Up + Op. figure 14. 

The percent of Np + 0g along each 
curve is constant (jahn) . 
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9 11 12 CH^ 

Figure 13.- Normal velocity of flame 

propagation in methane- 
air mixtures (Coward and Hartwell) . 




Figure 15. 
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Pigs. 16, 17, 18,19. 20, 2i 
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Figure 18 



Figure 19 
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Figure 20 
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^co+Og^ 00 m > 2 -Pco+o 2 =150 m > 
3 -Pco+o 2 =200 m ' 4 -Pco-j-o 2 =30 ° m > 

5-p C0+02 =760 MM 
Figure 21 



